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Sb- and Mn-activated calcium halophosphates were produced and their fluores- 
cence was investigated by the spectral method in its dependence on the activator con- 
centration and the weight ratio of the halogens. The trichromatic colour coordinates 
were computed from the experimental intensity distribution. The maximum quantum 
yield of the fluorescence as a function of the concentration of activators and halogens 
and a displacement of the fluorescence maximum of Mn was found to occur in the spectral 
range of 5760—5935 À when fluorine had been entirely replaced by chlorine. The flu- 
orescence mechanism is discussed in detail; the one consiting in direct resonance energy 
transfer from the sensitizer to the activator was found to agree best with the experimental 
results. According to the author's results, tbe range of the Sb to Mn 'energy transfer 
amounts to 50 cation loci. Perturbation theory as applied to resonance transfer of energy 
in solids yielded an explanation of the possibility of Mn transitions into forbidden states 
through the agency of the sensitizer, the return wherefrom gives rise to the orange flu- 
orescence characteristic of Mn. 


Introduction 


Calcium halophosphates were discovered by Mc Keag and Ranby in 1942 in Great 
Britain (1942). The basic substance of formula 3Cag (PO), Ca (FC), exhibits apatite 
structure; the activators are antimony and manganium. The compounds are intensely 
fluorescent when excited with the 2537 À mercury resonance line, and hence are 
widely applied in lighting techniques in low-pressure mercury lamps. They are of 
a high degtee of interest to scientific research, too, as belonging to the group of doubly 
activated luminophors, the discovery of which occurred at a rather late period. 

For this reason, the compounds have been the object of investigation by various 
authors. At first, principally the technology of production of calcium holophosphates 
was dealt with; later investigation centered on the mechanism of luminescence of the 
luminophors. 

Jenkins and co-workers (1949) elaborated a technological method of obtaining 
calcium halophosphates, and investigated the fluorescence of the latter in its depen- 
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dence on the concentration of the activators and on the atomic weight of the halogen 
used. 

Further, extensive work on calcium halophosphates was carried out by Butler 
and Jerome (1950). 

Bandel and co-workers (1951) investigated the properties of the fluorescence 
of calcium halophosphates as dependent on the chemical composition of the lumino- 
phor, the weight ratio of the halogens, and the amount of activators introduced. 

Doherty and Harrison (1954) worked out an extensive method of preparing 
the salts composing calcium halophosphates, and investigated the fluorescence of 
a number of luminophors obtained therefrom. 

It was the aim of the present paper to obtain a variety of calcium halophosphates, 
to investigate their fluorescence spectrum by the spectral method as dependent 
on the concentration of the activators and on the quantitative ratio of the halogens, 
to compare the results with those obtained previously, and to determine the mechanism 
of light emission. 


1. Technology of Production of Calcium Halophosphates 


The present author obtained calcium halophosphates by the method of sintering 
mixtures of the following salts: 


. CaHPO, 

CaCO, 

. (NH4, HPO, 

CaF, 

CaCl, 

Sob 

. Mn; (PO), > 


NAN PWN D 


These were obtained from the Production Centre of Chemical Reagents (Wy- 
twórnia Odezynniköw Chemicznych) at Gliwice, Poland, with the designation as being 
„of special purity for halophosphates“. Previous to sintering, the salts were throughly 
dried at about 300°C and ground to a grain size of 5 — 10 w. On mixing the compo- 
nents in the appropriate wieght ratio, the mixture was sintered during 35 minutes 
in an open quartz tube at 1050° to 1110°C. Mixtures containing less Mn required 
a lower sintering temperature, whereas those with a higher Mn content required 
a higher temperature. Previous to introducing the tube into the oven, the latter 
was heated to the temperature required. Sintering proceeded in an atmosphere of the 
gases leaving the tube, as: CO,, P,0;, NH,, N, and SbCl. The technology of the 
sintering process is known to consist on the following: the high temperature leads 
to dissociation of the different molecules of the mixture into ions, which then regroup 
to form crystals of calcium halophosphate with ions of the activators incorporated. 
The remaining ions combine to gas molecules and leave the site of the reaction, vola- 
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tilizing process, from 22 to 24% of the mixture volatilized. It is required that the 
atmosphere of sintering should not have an oxidating effect, as the activators can 
incorporate into crystal lattice of the calcium halophosphate in the form of ions 
of the lowest valences. This requirement is set by the ionic radii. The product, on 
cooling in an anhydrous atmosphere of very low oxygen content, was powdered in 
a mortar, thus yielding the ready luminophor. In this manner the author prepared 
two kinds of luminophors: one containing luminophors of the constant fluorine 
to chlorine weight ratio of 49 : 51 and a variable weight ratio of the activators, and 
another consisting of luminophors of a constant weight ratio of the activators and 
a variable ratio of the halogens. The foregoing amounts of activators and halogens 
in weight per cent relate to the mixtures previous to sintering. Since Sb and Cl, vola- 
tilize in the process of sintering in the form of SbCl, the respective figures serve 
primarily to provide orientation as to their amounts present in the luminophors.!) 


2. Account of Experimental Method and Apparature Used 


The fluorescence was investigated with an HCII—51 glass, three-prism spec- 
trograph. The various elements were assembled as shown in Fig. 1, wherein: 

Q — denotes the light source, 

D — the diaphragm limiting part of the light flow, 

F — Schott and chlorine filters, 

L,— a convergent quartz lenz, 

H — a container for holding the luminophor, 

L,— a convergent glass lens, 

S — the spectrograph slit. 


Fig. 1. Scheme of apparatus used in investigation. 


The luminophor T illuminated laterally with a TUV— 30 Phillips bactericide 
lamp. The exciting light vonsisted of the mercury line of 2537 Ä wavelength, the 


1) The true content of chlorine and the activators in the luminophors is at present the obiect 
of investigation.. ~ 3 
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remaining lines being absorbed in the Schott and chlorine filters. The spectrograph 
slit width amounted to 0.1 mm. The time of illumination ranged from 2 to 20 minutes. 
The lines obtained from a helium filled tube served for reference. Polish made ^UI- 
trapan Super“ plates were used for the photographs. The luminophors were inves- 
tigated at the constant temperature of — 20°C. 


3. Determination of the Intensity Distribution 


To this aim, two series of photographs were taken, one with a graduated atten- 
uator for transverse photometrisation, and one without the attenuator, for longitudi- 
nal photometrisation. On plotting the curves of characteristic blakening, the int- 
tensity distribution was determined for the wavelengths along which transverse pho- 
tometry of the plates had been carried out. Longitudinal photometrisation was used 
for interpolation. On one and the same plate were photographed: the lines of ref- 
erence, the spectrum of the normal lamp (which was the control spectrum), and the 
fluorescence spectrum of the respective group of luminophors. The plates were 
photometrized on a Kipp automatic microphotometer. The relative distribution of 
the intensiteis of the normal lamp spectrum was determined from Wien’s formula 


(Ornstein. 1958) 


1(A) = C, 475 exp | "T =) (3.1) 


wherein / (A) is the radiation intensity for a given wavelength A, and C,, C, are con- 
stants; 
C, does not intervene in comparative measurements of intesity ratios, 
C, = 1.438 cm. grad (a universal constant), 
T is the temperature of the colour in the Kelvin scale, and was determined as T = 
= (2840+10)°K. 

The method used in the present investigation for determining the intensities 
eliminated the spectral sensitivity of the plate. 


4. Determination of the Trichromatic Coordinates 


From the intensity distribution, the trichromatic coordinates x, y, z were deter- 
mined for the various luminophors; this —yielded an objective measure of the colouring 
of fluorecsence, which plays an important in lighting techniques. 

The following formulas were used: < 


x = f s (ZAI) da; y'= foam (2) ¥ (A) TA) das z — f Frum (2) 2 (A) 1A 
| (4.1) 
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wherein x(A), y(A) and z(A) are the coordinates of the isoenergetic source, and IA) — 
the intensity of source A of colour temperature 2854^K. 
X, y, z as functions of x’, y’, z' are given by the formulas 


, , , 


x y z 
AARON E EDO RC oe ET EUR . 
tay abe " Pe uz a a EST ODE C 
with 
x+y+z=l (4.3) 


5. Account of the Fluorescence in the Various Groups of Luminophors 


a. Group I. 


Group I contains 15 antimony activated luminophors, the amount of Sb increa- 
sing from 0 to 7.5%. Luminophor I-1 contains no activator and does not emit light 
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Fig. 2. Example of curves of fluorescence inensity distribution in luminophors of group I, versus the 


wavelength. Mn = 0%. I-2 Sb = 0.05 %, I-7 Sb = 1.5 %, I-15 Sb = 7.5 %. 
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on excitation. The distribution of fluorescence of the other luminophors over the 


various wavelengths is shown in Fig 2. 
The graphs show that the luminescence is blue, with maximum intensity at the 
wavelength of 4472 A. The dependence of the intensity maxima on the Sb concentra- 


tion is illustrated in Fig. 3. 
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Fig. 3. Maxima of fluorescence intensity versus Sb concentration Mn = 0 95. 
b. Group Il. 


Group II comprises 15 Sb- and Mn-activated luminophors, the manganium 
concentration being kept constant at 0.9 96, and that of antimony increasing from 0 
to 7.5 99. Luminophor Il-l contains no Sb and does not emit light on excitation 
with light of a wavelenth of 2537 À; however, on for ultraviolet excitation of about 
2000 À, it emits the orange fluorescence characteristic of Mn. The absorption spectrum 
of manganium activated caleium halophosphate lies in the spectral range beyond 
2200 À, and thus radiation of higher energy, as e. g. X-rays, cathode rays and far 
ultraviolet can be used for excitation (Botden and Krüger 1948). Luminophors of this 
group exhibit two ranges of fluorescence: the orange one of Mn, and the blue due 
` to Sb. About 5340 À, the emission band of the two activators overlap. The intensity 
versus the wavelength is given in Fig. 4. 
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Fig. 4. Example of curves of fluorescence intensity distribution in luminophors of group II, versus the 
wavelength. Mn = 0.9%. II—2 Sb = 0.05 %. II—8 Sb = 1.88 %. 


The distribution of the Sb and Mn fluorescence intensity maxima, versus the rate 


of Sb introduced, is shown in Fig. 5. 


c. Group II. 


The third group comprises 13 antimony- and manganium-activated luminophors 
containing the constant amount of 1.88 96 Sb, and wherein Mn varied from 0 to 4.05 95. 
The intensity distribution versus the wavelength is shown in Fig. 6. 

The fluorescence consists of three bands: one deriving from the manganium 
and having its maximum in the same spectral range as group II, and two belonging 
to the antimony exhibiting maxima at 5016 À and 4336 À. The manganium emission 
intensity increases steeply with the amount of this activator, and attains its maximum 
value in luminophor III-8 at 1.8 % Mn, whereafter it decreases slowly. The antimony . 
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Fig. 5. I Graph of fluorescence intensity maxima of Sb versus the Sb concentration. Mn = 0.9 %. 
II. Graph of fluorescence intensity maxima of Sb versus the Sb concentration. 


bands are of a different charcter: radiation of longer wavelength appears first in lumi- 
nophor III-3 and persists as far as luminophor III-10. Its intensity decreases 
rapidly. The band of shorter wavelength appears in all luminophors. The decrease 
of the fluorescence intensity occurs exponentially. The distribution of the intensity 
maxima, as dependent on the amount of Mn incorporated, is shown in Fig. 7a and 7b. 


d. Group V. 


The fifth group contains five luminophors of the constant activator contents of: 
Sb = 2.25 %, Mn = 0.9 9$, and of variable fluorine to chlorine weight ratio. 

The intensity distribution versus the wavelength is shown in Fig. 8. 

The distribution of the intensity maxima of the Mn and Sb emission, as versus 
the weight of the halogens, is shown in Fig. 9a and 9b. 

From Fig. 8, a displacement of the Mn fluorescence maxima in the range of 
5764—5935 A is clearly seen to occur when fluorine is entirely replaced by chlorine. 
However, this causes the maximum of Sb emission to displace but insignificantly 
towards longer wavelengths. The strongest intensity of Mn and Sb fluorescence occurs 
in luminophor V-2 of fluorine to chlorine weight ratio 2 : 1. 
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Fig. 6. Example of curves of fluorescence intensity distribution in luminophors of group III, versus 
the wavelength. 


IIlI—-4 Sb = 1.88% III-8 Sb= 1.88 9$ III—11 Sb- 1.8896 
Mn = 0.6 95, Mn = 1.8 95, Mn = 3 9$. 


6.Theoretical Interpretation of Luminescence Mechanism 
. in Luminophors Investigated 


A crystal wherein transition of an atom or ion into an excited state has occurred 
‘by absorption of radiation or otherwise, can return to the normal state by ridding 
itself of the energy absorbed. The latter usually is dissipated in the form of heat. 
In some cases, however, it is re-emitted in that of radiation, and the substance ex- 
hibits luminescent properties (Mott 1956). It is seldom that these are exhibited by pure 
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substances. This is the case of manganium halogenides, rare earths salts, platinocya- 
nides, wolframates, molybdates, and uranyl salts. Randall (1939) suggested that 
a pure substance is fluorescent if each elementary cell of the crystal contains an ion 
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hole should be shielded off from the lattice vibrations. It is here that the crystal lattice 
defects intervene. At the usual temperatures no ideal lattices exist, since defects 
occur therein as the vibrations can cause an atom or ion of the crystal to depart from 
its usual lattice point and to migrate to an interstitial site or to the crystal surface. In 
certam circumstances, foreign atoms or ions, termed activators, can be incorporated 
at the void crystal points, forming centres wherein the electron and hole recombine 
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Fig. 7b. Fluorescence intensity maxima of Mn, as versus concentration of Mn. Sb = 1.88 95. 
with the emission of radiation. The reason for this is that, in a manner hitherto not 
fully elucidated, the electron and hole are better shielded off at the conters than else- 
where from the effect of the lattice vibrations and thus the probability for them to 
transfer their energy to the surrounding atoms or ions in the from of heat is lesser. 

The conductivity band is essential for the mechanism of luminescence of lumino- 
phors having a single activator. Two cases can occur here (Jabłoński 1950): 

1. Absorption of exciting radiation is related to the transition of an electron 
from the activator level to that of conductivity. Emission occurs when one of the 
electrons of the conductivity level (an excited electron) falls to onto the activator 
level (Fig. 10). 

2. Excitation occurs by the transfer of an electron from the invested fundamental 
band to the conductivity band. An electron from the invested activator level occupies 
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Fig. 9. II. Fluorescence intensity maxima of Sb, as versus the concentrations of ‘fluorine and chlorine. 
II. Fluorescence intensity maxima of Mn, as versus the concentrations of fluorine and chorine. 


primary activator or sensitizer, and manganium the secondary activator or activator 
proper. Such luminophors emit what is known as sensitized fluorescence. According 
to the general theory of sensitized luminescence, the activators replace cations of the 
basic substances (Williams 1954). Thus, the ions of the two activators replace calcium 
cations. Manganium and antimony are incorporated in the lattice in the form of 
bivalent and trivalent ions, respectively. Together with the antimony ion, an oxygen 
ion is incorporated replacing a halogen ion. : 
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In the present investigation, the luminophors were excited with the mercury 
9537 A line. The role of antimony as sensitizer in this group of luminophors consists 
in bringing the absorption spectrum into the range of excitation, since neither the 
basic substance by itself nor when manganium-activated absorbs within this range. 
Part of the energy absorbed by the sensitizer is emitted in the form of its specific 
fluorescence, while part is transferred to the activator and emitted in that of the 
emission spectrum characteristic of the latter, i. e. of manganium. The fluorescence of the 


Fig. 10 l Fig. 11 


Fig. 10. Diagram showing fluorescence mechanism in electron transition from activator level to con- 
ductivity band. 
Fig. 11. Diagram showing fluorescence mechanism in electron transition from fully invested band to 
conductivity band. 


luminophor is always complex, and consists of the emission bands of either activator, 
the ratio of the intensities depending on the concentrations of the latter. 

A problem of far reaching significance arises at this stage, namely, that of the 
transfer mechanism of energy from the sensitizer to the activator. The literature 
brings a number of theories aiming at an explanation of the phenomenon of sensi- 
tisation. Some had to be rejected as being at variance with the experimental facts. 
Three examples will be given. 

l. The hypothesis of energy transfer from sensitizer to activator through the 
mediation of the conductivity band had to be rejected, since experiments proved 
that sensitized luminophors do not exhibit photoconductivity. 

2. The hypothesis of energy transfer from excited sensitizer ion to activator 
ion over the ground state or excited state of the crystal lattice of the principal suk- | 
stance had to be rejected, too, as requiring a quadratic relation between the quantum _ 
yield of fluorescence and the intensity of the exciting radiation. 

3. The hypothesis of the formation of complex sensitizer and activator centres, 
termed coactivators, is also at variance with the experimental results (Botden 1952). 

~The most plansible theory consists in direct energy transfer from sensitizer to acti- 
vator in a resonance process resembling the one occurring in liquids and gases. We | 
now introduce the following two assumptions: - | " 
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1. A random distribution of ions of the sensitizer and activator throughout 
the luminophor is assumed, and 

2. Resonance energy transfer between sensitizer-sensitizer and sensitizer-acti- 
yator ions is assumed (19522). Energy is transfered from an excited ion of the sensitizer 
to an activator ion not only by direct resonace involving the two ions, but also over 
other ions of the sensitizer. Hence, energy absorbed by the sensitizer at a considerable 
distance from the activator can be transferred to the latter by other sensitizer ions 
(“enrgy transfer“). 

We shall now proceed to derive a formula for the energy transfer efficiency. 
The latter is defined as the ratio of the quantum fluorescence yield of the activator 
and the sum of those of the sensitizer and activator: 


A 
Ze ER (6.1) 

wherein Z denotes the energy transfer officiency, 

S — the quantum fluorescence yield of the sensitizer, 

A — the quantum fluorescence yield of the activator. 

Let us consider all possible paths of energy transfer from sensitizer to activator, 
namely: 
a) an S ion that can transfer energy directly to A (the Sy — A process). 
b) an S ion whence energy can be transferred to A exclusively RUE the inter- 
vention of another S ion (the S, > 5, > A process). 
c) an S ion from which energy transfer to A can occur only over two S ions (the 
Sp 2 Si — Op — A process). 

Let k denote the number of cation loci surrounding 5 u having the property 
that energy can be transferred from S to 4 if 4 occupies one of the k-th loci. 

The probability for energy transfer is given by the circumstance that whereas 
absorption of excitation energy at S situated at any one of the k loci from A involves 
fluorescence of A, absorption at an S further remote can lead to fluorescence, of S. 

Let X, and X, denote the ratio of the number of ions of the activator and sensi- 
tizer, respectively, to the number of cation loci in the crystal. By assumption (1). 
the probability for finding A at a given distance from S is X,, whereas the probability 
that A is not situated at this particular point is 1—X,. Then P, = (1—X,)* means 
that at each of the & loci not A but some other ion is situated, and P, = 1-(1—X,)* 
is the probability for finding A at the k-th locus. 

The contribution to Z from the Sy — A process described in (a) is in by the 


ratio Pa: (P, + Pa), 1. e. by i 
1—(1— X,)* (6.2) 


Let Z denote the number of cation loci surrounding 5 and having the property 
of energy transfer between one S and another S situated at one of the /-th loci, and 
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p — the probability of such transfer. Then the contribution to Z from the 5, > S> A 
process of (b) is given by the expression. 


(cL A A (6.3) 
I Ji III 


wherein I is the probability that there is no A at any of the k-th loci surrounding the 
initially excited Sọ II is the probability for S, to be present at one of the k-th loci 
around Sọ and for energy transfer from Sp to S,, and III — the probability that A 
is to be found at the k’-th locus in the neighbourhood of S,, and that the energy 
received by S, is transferred to A. Expression III is analogous to that of (6.2) wherein 
k has been replaced by k’, since some of the k’ loci around S, can be among the k 
loci surrounding Sọ; the number of overlapping loci is (k—k’). 

Similarly, the contribution to Z of the Sy > Sı > S > A process (c) is given by 


(1X) p [I (1—3X)10— x)" pll— E E X eto 
I | II III 


wherein I is the probability for Sọ — S} energy transfer from a given locus / in the 
neighbourhood of S, and for no A to be present at any of the loci k or k’ around S, 
or Sy II is that for the presence of an S, at one of the /-th loci around S, and for the 
occurrence of 5$, — S, energy transfer, and III — the probability for an A to be situ- 
ated at one of the k”-th loci around ie accounting for overlapping of loci of this type 
and those of types k and k’ around S, and $,. 

Similarly the contribution to z from G S PS ass ‘Ss — A process 
etc., can be computed. 


Adding expressions (6.2), (6.3) and (6.4), we have 


k(k— 1) 
kX,— —5 


; k(k— 1) (k —2 es 
Xi + plk po ue M [Exe +) = 


KE 1 
ran! 24 pik’ | Ar (6.5) 


If. X, is small with respect to 1/k and X, is small with respect to 1/], the terms 
in (6.4) decrease with higher powers of X, and X,. 


Denoting by Z,, Z and Z, consecutive approximations to Z, the index equalling 
the maximum value of the sum of the powers of X, and X,, we have 
Z = kX, (6.6) 
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By eq. (6.6), the quantity Z does not depend on the concentration X, of the 
sensitizer, and is constant for a given concentration X, of the activator. Hence, 

Z 

k =a EE 


T (6.9) 


This state can be attained if the concentration of the sensitizer is high as com- 
pared to that of the activator and all the activator ions lie within the reach of the 
sensitizer. It is then that direct Sọ — A energy transfer takes place. If this can be 
attained experimentally, thequantity k can be computed from eq. (6.9). 

By eq. (6.7), Z is seen to increase linearly with X, if X, is kept at a constant 
value. Substituting X, — 0 in eq. (6.7), we have 


k(k— 1) 


2 
"NE 


Z(x,-0) = kA, 


whence 
_ 1—[(l— 22) xL g]* 
— X; 


From eq. (6.10), k can be computed. In practice this is possible if the concentration 
of the sensitizer is small as compared with that of the activator, and not all the ions 
of the activator are within the range of action of the sensitizer ion. Then energy transfer 
takes place through the agency of another sensitizer ion, according to the process 
EN a A, i 

By eq. (6.8), Z is a quadratic function of X,. This corresponds to a process 
wherein energy can be transferred from S to A exclusively over two S ions, i. e. accor- 
ding to Sy > Sı > Sa > A. This should be expected to occur at the very lowest 
concentrations of the sensitizer. 

The foreging theoretical considerations will now be applied to calcium halo- 
phosphate luminophors. 

. By simple considerations, the efficiency of energy transfer, Z, equals the ratio 

of the fluorescence intensity of the activator and the sum of those of the activator, 


k (6.10) 


and sensitizer: 


IMn 

RER DU 
with Jmn and Ig, denoting the intensities of emission of Mn and Sb, respectively. 
It result that, if the distribution curves of the intensities are known the quantity Z 
can be computed. The numerical values, as dependent on the Sb concentration, 

are shown in Fig. 12. 
From Fig. 12 the curve of efficiency of energy transfer, Z, as dependent on the 
concentration of antimony, si seen to consist of three sections: section J, which cor- 
responds to the lowest Sb concentrations, is an ascending curve and accounts for 


^ 
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Fig. 12. Energy transfer efficiency versus Sb concentration. Mn = 0.9 %. 


the Sy — Sı > S, > A process of energy transfer; section II, for Sb concentrations 
ranging from 0.38 to 1.13 96, is an ascending straight line accounting for the 
Sy > Sı > A process; section III, for Sb concentrations of from 1.13 to 7.5 96, is a 
straight line parallel to the Sb concentration axis and corresponds to the Sy > A 
process of energy transfer. Thus it is to be concluded that the foregoing theore- 
tical considerations are in agreement with the experimental results. 

The value of k as computed by eq. (6.10) from section II of the curve amounts 
to 61, whereas that computed by eq. (6.9) from section III is 50. The latter is the 
move reliable of the values, seeing that former is obtained by interpolation of the curve 
down to X, = 0 and results from no more than 4 points. Botden (1952) obtained 
values of k for Mn ranging from 26 to 57, as computed from eq. (6.10), with several 
sensitized systems. However, Botden had insufficient experimental material at his 
disposal, and in none of the systems investigated did he attain independence of the 
efficiency of energy transfer with respect to the sensitizer concentration according 
to eq. (6.6). 

Dexter (1953), from theoretical computations of dipole radiation, obtained 
a value of amounting to 30 for Mn. 

Similar considerations hold for the III group of the luminophors containing 
a constant concentration of Sb and variable Mn. 

By eq. (6.6), Z should now be expected to increase linearly with the activator 
concentration. Results as computed for Z versus X, are shown in Fig. 13. 

The values of the efficiency of energy transfer as computed for the III group 
of luminophors corroborates the foregoing conclusions. For the Sb concentration, 
which is 1,88 96, energy transfer is seen to occur according to the equation Z = kX,, 
i. e. the So > A process. For Mn concentrations exceeding 3 96, transfer efficiency 
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Fig. 13. Energy transfer efficiency versus Mn concentration. Sb = 1.88 %. 


decreases. Presumably what is known as fluorescence quenching due to activator 
concentration occurs heve, i. e. the effective number of emitting centres diminishes. 
According to Schulman (1954), the activator can absorb but cannot emit when within 
the range of action of another activator ion. Since the fraction of activator ions within 
the range of another activator ion increases with the concentration, the fluorescence 
yield must necessarily decrease towards higher concentrations. 


7. Theory of Resonance Energy Transfer in Solids 


The significance of sensitized luminescence in inorganic solids resides in the 
fact that resonance energy transfer can take place between a permitted transition 
in the sensitizer and a forbidden transition in the activator. Thus the activator in the 
present instance manganium yields luminescence subsequent to the transfer of the 
energy absorbed by the sensitizer, since, by itself, it does not absorb directly. Accord- 
ing to Dexter (1953), the entire process of resonant transfer of energy comprises 
the following five stages: 

l. absorption of a photon of energy E, by the sensitizer, 

2. relaxation of the crystal lattice around the sensitizer ion, to such an extent that 
the available electron energy in radiative transfer from the sensitizer shall be E, — Ey, 

3. transfer of the’energy P, to the activator, 
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4. relaxation of the crystal lattice around the activator ion to such an extent 
that the available energy emitted in radiation shall be Ej < £i, 

5. emission of energy P. 

Stages 3 and 4 require 10-13 sec each, or maybe longer, for excess energy to be 
dissipated by the phonons over the crystal lattice, leading to Stokes’ law. Stage 3 takes 
a time that depends on the distance of the surrounding activator ions, i. e. on the 
concentration X,. Stage 5 requires 1078 sec or much more, according to the degree 
of forbiddenness of the transition. Since Stokes’ law is usually fulfilled, the activator 
represents a kind of “trap“ for the electron excitation energy, i. e. the excited activator 
cannot enter into resonance with the sensitizer (E; < Ej) and no energy transfer 
can take place from activator to sensitizer. The foregoing considerations serve to 
explain two facts of the domain of sensitized luminescence, namely: 

1. in all sensitized systems, fluorescence of the activator exhibits a longer wave- 
length than that of the sensitizer; 

2. the fluorescence decay time of the sensitizer is much shorter than that of the 
activator, e. g. amounting to 10-9 sec for Sb and to 107? sec for Mn. 

We shall now consider resonance energy transfer from a theoretical point of view, 
starting from the perturbation theory (Herzberg 1939) of quantum mechanics. By 
this theory the magnitude of the pertubation between two states depends on the 
matrix element W,; and on the differences in energy between the non-perturbed 
levels. The matrix element ]/,; can divided in two parts: one dependent on the 
electron and rotatory parts of the wave function, and one dependent on the oscil- 
latory part thereof. The latter part is 


W,— f| ye W yt" dr (7.1) 


wherein y% and y?" are oscillatory wave functions in Schrödinger’s equation, and W” 
is a function of the distances between the nuclei. By eq. (7.1), a strong perturbation 
occurs only when the oscillatory parts of the wave functions overlap appropriately. 
This can occur only for two oscillatory levels situated in the neighbourhood of the point 
of intersection of the potential curves. In luminescennce, a perturbation between 
two oscillatory states appertaining to two different ions (or atoms) leads to energy 
transfer from one of the ions (or atoms) to the other. In the perturbating state apper- 
tains to an ion of the impurity, usual luminescence quenching occurs. In sensitized 
luminescence, an ion of the activator proper is the perturbating factor. At the mo- 
‘ment of excitation it remains in the normal state, as transition to an excited state is 
forbidden by the selection rules. However, if it possesses excited states of energy 
almost equal to that of the excited states of the sensitizer, the energy of the latter 
is transferred to the activator, raising it to an excited state. Hence, in sensitized 
luminescence, the roles of sensitizer and activator proper can be played only by ions 
fulfilling the two conditions: 


l. they must exhibit excited states of almost equal energy, 
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2. for obtaining maximum efficiency of energy transfer, the differences between 
the various excited energy states, both in the sensitizer and activator, should 
be small. 


Numerous tests proved that these are best fulfilled by the combination of Sb 
as sensitizer and Mn as activator proper. Antimony can be replaced by arsenic, 
bismuth, lead or tin, thoung the result is so satisfactory. The only activator proper 
known in sensitized luminescence is, however, manganium. 

Applying the perturbation theory outlined above to sensitized luminescence 
in solids, the following conclusions may be deduced: 

l. Two oscillatory states of sensitizer and activator, of almost equal energy, 
interact in such a way that energy from the sensitizer is transferred to the activator. 


r 


Fig. 14. Curves of potentional energy of fundamental and excited states for antimony and manganium. 


2. Maximum interaction occurs if the potential energy curves of either. state 
approach one another or intersect. 

3. Hence results the relation between the potential energy curves in the fun- - 
damental and excited states for Sb and Mn shown in Fig. 14. 


U —- potential energy, 


r —— distance between nuclei, 

N — fundamental state, 

Pyn— excited state in Mn, 

Ps, — excited state in more remote Sb, 
P$, — excited state in nearer Sb. 


Energy absorbed by the Sb ion very remote from the Mn ion is transferred 
to a nearer ion, whence it is, in turn, transferred to the Mn ion. 
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ON THE POSSIBILITY OF DETERMINING THE DENSITY OF 
DISLOCATIONS BY MEANS OF AN X-RAY SPECTROGRAPH 
WITH OSCILLATING FILM*) 


By Teresa BEDYÁSKA 
Institute of Physics, Polish Academy of Sciences, Warsaw 
(Received Nowember 17, 1959) 


A method of computing the distribution fanction of mosaic blocks im the ars 
Í is proposed. For this, it ie required that the forme of Bragg’s reflex, s obtained on am 
l X-ray spectrograph with fixed film, and that of the reflex, as obtained on 2 spectrograph 
with oscillating film, be known. If the distribution function of mosaic blocks ia known, 

the density of dislocations can be computed. 


1. Introduction i 


The boundaries between blocks in a erystal of mosaic structure are formed by 

. the dislocation lines. If the crystal consist structurally of blocks of edge length t 
and if the distance between the dislocation lines on the boundary separating the 
blocks is h, the density of dislocations on the boundary proma — !ht? = 1/ht (Gar, 
Hirsh, Kelly 1953). If the blocks in the crystal were formed by a system of parallel 
dislocation lines having one and the same Burgers’ vector, the density oramg t computed 
from Gay, Hirsh and de ee ro gtd would be numerically equal to the nunäber 


cuti 
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micro-beam technique, reflexes from the various blocks can not be obserwed separatly, 
Gay, Hirsh and Kelly (1953) assume a system of blocks with a single dislocation 
per one inter-block boundary. With this assumption, the total dislocation density is 


given by formula: 
2 
a 
oc (=) (2) 


In order to find the dislocation density, the mean angle between the blocks 
should be known. The latter may be determined if the angular distribution function 
for the mosaic blocks in the crystal is known. This function can be obtained from 
analysis of the line form yielded by the X-ray spectrograph with oscillating film 
(Auleytner 1958, 1959 a, b), and comparing it with the one yielded by a spectrograph 
with immobile film. 


2. Principles of Method of Determining Mosaic Blocks Distribution Function | 


The principle of the method is given in Fig. 1. An ideally parallel and narrow 
beam of X-rays is incident at point O on the ideally smooth crystal surface. The 
rotation axis of the crystal is at point O. The crystal rotates about axis O ang displaces 
parallel to the reflecting surface. On the immobile film F, a line of the form given 


l 
| 
| 


On the Possibility of Determining the Density 445 


In experimental investigation, the effects of: 1) divergence of the beam, 2) width 
of the slit, and 3) inequalities of the crystal surface and displacement of the axis 
of rotation with respect to the reflecting surface should be accounted for. 


3. Effect of Divergence of Beam in Horizontal Plane on Line Forms Obtained 
on Immobile and Oscillating Films 


In the case under consideration, a beam, divergent in the horizontal plane (the 
plane parallel to the rotation axis of the crystal) and parallel to the vertical plane, 
is incident on the crystal K from an ideally narrow slit S (Fig. 2) on the focussing 
circumference. The form of the divergent beam is described by the function (Rf) 


Fig. 2. Effect of idivergence of the beam in horizontal plane on forms of the lines. 


with f denoting the angle subtended by the beam SP incident on point P of the 
crystal and the direction of the beam SO incident on its rotation axis, and R denoting 
the distance of the slit S from the rotation axis O. Let œ denote the desorientation 
angle of the mosaic blocks. The displacement, on the immobile film, of the beam 
reflected at point L (Fig. 2) with respect to that reflected at O, as resulting from the 
mosaic, is a: . 

r a = 2 cos Ÿ t sin a 
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wherein @ is the angle fulfilling Bragg's equation, and t — the distance of the refllec- 
ting point L of the crystal from the axis O: 


R 


Se 


whence 


cos 4) 
sin ( + æ) 
As the crystal K rotates, the angle between its surface and the direction SO varies. 
Let o be the angle by which the one subtended by the direction SO and the crystal 
surface differs from 0. The condition for the fulfillment of Bragg's equation for a block . 
deflected by æ is given by 


a= sin f sing — 


The line form, for immobile film, and rotated and displaced crystal, is given by the 
function G(x): 


-Foo -Foo N 
1 cd , TEDER 
cw- f fe (on Rinne a) x 
x ®(RB) {RP — Rf) d (Rp) d (RB) Tous 
The line form, for mobile film, and rotated and Dre crystal is given by the 
function M(x): D 
T ; 
os 9 
H (x) = WEL EN 
(x) / IE [-- 2R tm í 7 E sin fs sin (p — f) — Re|x 
-00 —00 


uni zn alti vn 
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wherein 


h(x) = ji g (x — Ra) f (Ra) d (Ra) 


The function A(x) is the one that would result on the oscillating film with a parallel 


beam eq. (3) 


4. Effect of Finite Dimensions of Slit on Line Forms G(x) and H(x) 


If the slit of finite. width, the repartition of the, intensity thereon is accounted 
for by the function Ø, (z), From each point of the slit S (Fig. 2), a beam diverges. 
The function Ø, (z-- RB) accounts for the divergence in the horizontal plane of the 
beam leaving point z of the slit. With such a slit, the line form on the immobile film 


becomes 
Loo 
G(x) — f g (4 (æ — 2)) D, (2) dz (8) 
and, on the oscillating film: 
+00 +o 
H(x)= f f h(- RB—2 D, (2) D, (€ + RB) dz a(RP) 
+00 ta - 
= if f h (x — t) ®, (t) ®, (2) dz dt 
whence 
H(x)= f h(æ— i) D, (t) dt (9) 
wherein 


+00 
D, (t) = f ®, (2) 9, (0) de 


The function @, (t) describes the form of the beam leaving the slit, at a distance there- 
from equal to that of the crystal from the slit. 


5. Effect of Inequalities in Crystal Surface and of Displacement of Rotation Axis 
with Respect to Crystal Surface 


Inequalities inherent in the crystal surface, and displacement of the axis of rot- 
ation with respect to the crystal surface can result in some additional widening of the 


lines G(x) and H(x). , 
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Displacement of the crystal surface by the amount t causes the line g(x) on the 
film to displace by d (Fig. 3). g(x) describes the line form that would arise on tho 
immobile film with an ideally narrow and parallel beam: 


d = 2 cos 0 £ (10) 


Displacement of the rotation axis of the crystal with respect to its surface by the 
amount / gives rise to the variation tọ in the distance of the reflecting surface during 
rotation of the crystal (Fig. 4). If y denotes the angular width of the line reflected, 
then 


tery cp T as ole C (11) 


Such difluence in G(x) and H(x), arising from this variation in the distance of the 
crystal surface, can be neglected. 


ft 


hiia cie did 


—"-———— 
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If the widening of the line on the immobile film, as resulting from inequalities 
in the surface scanned by the X-ray beam when the crystal is not displaced, be neglec- 
ted, the widening resulting from surface inequalities over a larger area, as investi- 
gated when the crystal is subjected to displacement, can be eliminated. In order 
to eliminate this latter widneing, it is necessary that, in addition to the line forms 
f(x) obtained with displaced crystal and oscillating film, as well as G(x) obtained 
with displaced crystal and immobile film, the line form obtained with non-displaced 
crystal and immobile film should be known. If the surface exhibits inequalities, the 
line resulting on the immobile film when the crystal is subjected to displacement 
assumes the form 

+00 
G(x) = f XO C @—Hae (12) 
— oo 
The one obtained on the oscillating film as the crystal is subjected to displacement 
assumes the form 
+00 
H(x)= f X() B'(— i) dt (13) 
—oo 
The functions G'(x) and H'(x) describe the line forms that would arise from the 
same crystal, assuming its surface to be smooth. Since the inequalities within the limits 
of the beam, at non-displaced crystal, are neglected here, the function G'(x) will 
be identified with the function obtained when the crystal is not subjected to displace- 
ment. As the functions G'(x), G(x) and H(x) are known experimentally, the unknown 
function X(x) of the widening arising from surface inequalities can be eliminated. To 
determine the function of angular distribution of the mosaic blocks, use is made of 


eqs. (3), (8), (9), (12) and (13): 


+00 
h(x) = f fo) say) dy 


Too 
GG) = f && «—»» 9,0) dy 


+00 
G() = [ XW) G'&—» dy 
+00 


H' (x) = f h(x—y) 9,0) dy 


0X 


+00 
HQ) = f X0) H'@—y) dy 


"y 
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The foregoing system can be replaced by the single equation 


+00 
H* (a) = f censo (14) 
together with the four integrals 

+00 

H* (x) f Hx — y) F, yy EN (15) 
+00 

G(x) = f Ce—y) Fy) dy (16) 
-4-00 

EQ) = | 6 @—y) D, (27) dy (17) 

+00 
FQ) = f € G&—»)9,0) dy (18) 


f(y) is the distribution function of mosaic blocks required. The following function- 
are given experimentally: H(x), G(x) G'(x), B(x) and o, (x) 

H(x) accounts for the line form, at displaced crystal, on oscillating film. 
G(x) accounts for the line form, at displaced crystal, on immobile film. 

G'(x) accounts for the line form, at non-displaced crystal. on immobile film. 
D,(x) describes the form of the beam leaving the slit, on the latter. — 
B, (x) describes the form of beam leaving the slit, at distance therefrom equal to that 
of the crystal from the slit. = 


6. Te. er Divergence within Vertical Plane 
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wy _ | 
cos? ^ 
sin Ù cos 9 is ^y t) te y «t uy 
a 2 cos®# ' 8 cost} ' 16 cof ð ^" 
l tg?^g l tg*g lo iot. | 
1 — cos? 9 | — oes 
= | 2 cos? A. 8 costd ^ 16 costo =) 


at small values of 9 
tg x ~ $tg 0 tg? p (1 4- 31tg? p 4- 11g* p ...) 


a ~ gtg Ütg?g (20) 


If the beam exhibits divergence both in the vertical and horizontal planes, a line 
will result on the photographic plate, of a width determined by the horizontal and 
vertical divergences of the beam. 

Divergence can be removed by reflecting the beam from two immobile crystals. 

Using a point focus X-ray tube, it is possible to operate with a beam divergent 
in the vertical plane and to employ eq. (14) and integrals (15), (16), (17) and (18) 
in computing the distribution functions of the mosaic blocks. A point focus yields, 
at each point of the film, an image resulting from the beams exhibiting one and the 


Fig. 5. Effect of divergence of the beam in vertical plane of forms on the lines. 


same value of the angle p. The beam k leaving the focus at an angle @ in the plane 

perpendicular with respect to the beam normal to the slit and rotation axis of the 

crystal, fulfills Bragg's conditions if the angle between the crystal surface and ky 
amounts to 


§—a = 0 —1igóbtgig (21) 
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Variations of the angle a = % tg 0 tg? y with ® can be neglected within the width 
of the line investigated. These decrease as the angle 9 diminishes. Hence, for any 
beam k subtending a small angle with Kp, the line forms: G(x) resulting on an immo- 
bile film, and H(x) — as obtained on an oscillating film — undergo no modification. 
The lines G(x) and H(x) are displaced. In experimental work, the height of the tube 
focus and photometer spot should be made so small as to provide for divergence 
Ag in the vertical plane within the height of the photometer slit yielding much smaller 
displacement of the line than the dispersion of the mosaic blocks investigated. If 
tg p cos? pAg is considerably smaller than the angle of dispersion of the mosaic 
blocks, then divergence within the vertical plane will not have contributed towards 
deforming the lines G(x) and H(x). In this case, the dispersion function can be compu- 
ted from eqs. (14), (15), (16), (17) and (18). 


7. Effect of Widening from Reflexion at Interior Faces, and from Surface 


Inequalities 


Eqs. (14), (15), (16), (17) and (18) do not serve the purpose of eliminating wid- 
ening as arising from the surface inequalities present within the width of the beam, 
and from reflexion at crystal faces situated at a certain depth beneath the outer surface. 
This latter circumstance can come to play an important part in crystals exhibiting 
a low rate of absorption, but can be eliminated by measuring the Bragg line form 
on an immobile film disposed at one half of the distance from the crystal in mobile 
film investigation. The block distribution function is now computed from the 
following formulas: 


+00 
ha) — f fo)e(x—» dy 


+00 


G(x) = f g(x—y) 9,0) dy 
ds 
H(x) = f h(x—5) 9,0) dy 


With the foregoing system of equations, error arising through surface inequal- 


ities and reflexion at interior crystal faces is eliminated. The system of equations can 
be replaced by the equation 


Too N 
H*(x) = f C(x) f(y) dy © (22) 


together with the integrals 


+00 
H’ (x) = | Ho) ®,(x—y) dy | Q3 
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and 
+00 


G'(x) = | G(y) 9, (x — y) dy (24) 
f(x) is the distribution function of the blocks required. 
H(x) accounts for the line form on the oscillating film, at distance R from the crystal. 
G(x) accounts for the line on the immobile film, at distance 14 R from the crystal. 
D (x) describes the shape of the beam incident on the crystal, at distance 14 R there- 
from. 
D, (x) describes the shape of the beam incident on the crystal, at the latter. 


8. On the Form the Block Distribution Function as Resulting from Variously 
Oriented Systems of Dislocation Lines 


The presence of dislocations gives rise to a system of blocks within the crystal. 
The blocks formed by edge dislocations are mutually rotated about the dislocation 
line. Those formed by screw dislocations are rotated about the normal to the plane 
of sliding. 


Assume the rotation axis of the crystal to coincide with the y-axis of the system 
of reference (Fig, 6), the reflecting plane of the crystal with the xy-plane, and the 
incident beam to lie in the xz-plane. Then a beam S reflected from a block of surface 
perpendicular to the axis also lies in the xz-plane. A beam §, reflected from a block 
rotated by œ about the y-axis lies in the xz-plane and yields, on the mobile film, a reflex 
removed by the angle œ from that yielded by s. For a block rotated by y about the 
x-axis, the components of the unit vector of the incident beam Sp. and reflected beam. 
Sp are: 


ê Soa (— cos (8 + g), 0, sin (8 + g)) 
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with @ denoting the angle fulfilling Bragg's equation, 


sin? cos sind n ae 


sno cue x 
2 cos y cos y cos? d 


for small values of the angle y 


sin p = itg By? 


OW EHEGNMA r s2 
(ent END sin 9 272) 
cos y cos y 


The beam s, yields a reflex that is displaced by the angle Ad with respect to that 
of s in the direction perpendicular to the line on the film (in the xz-plane): 


sin 3? cos 2y a 


sin (9 + A9) = 4 sin? 9 sin? y)-% 


cos y 
hence, for small values of y 
. Ad ~ (2 sin? 9 — $) tg By? (25) 
The displacement along the line on the film (in the y-axis) is 


y' = sin 9 sin 2 y (cos? 0 — sin? 9 + sin? 0 sin? 2 y)1 
at small values of y 
y' — 21g 0y (26) 


By eq. (21), with a point focus X-ray tube, points of the line distant along the 
latter by 2 tg by are separated in the direction perpendicular to the line by the amount 
2 tg°9y?. Finally, the widening of the line due to the presence of blocks rotated by the 
angle y about the x-axis amounts to 


ô= (21g* 0 + 2 tg 0 sin? 0 — $ tg 9) y? (21) 


If y is of the same order of magnitude as the angle of maximum dispersion CM 
of the blocks rotated about the y-axis, then 6 is considerably smaller than € m: Hence, 
the deformation of the line due to blocks rotated about the x-axis can be neglected. 


A crystal block rotated by an angle g about an axis lying in the reflecting xy-plane 
of the crystal and subtending the angle g the with the y-axis (Fig 7) reflects equi- 
valently to one rotated by æ about the y-axis and by y about the x-axis. The angles 
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V, 9, & and y are related by the following formulas: 
sin y = sin y sin p (28) 
sin @ = sin y cos g (1 — sin? y sin? p)” (29) 
for small values of the angle y 
y=ypsn gy (30) 
a = ycos p (31) 


A block rotated by an angle y about an axis subtending an angle y with the 
reflecting plane yields a reflex upon the film such as would result from a block rotated 


Fig. 7. 


by the angle v, about the projection of the axis on the reflecting plane. The angles 
y, x and y, (Fig. 8) are related by 


in # — in © 32 
sinn = cos y sin -5 (32) 
at small values of y 
VET COP 2, 


At small y a block rotated about an axis subtending the angle x with the reflecting 
plane yields, on the mobile film, a reflex as would result from one rotated by æ about 
the y-axis: À 


a= y, cos p = y cos p cos y = pcos ® (33) 
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g denotes the angle subtended by the projection of the rotation axis of block on the 
reflecting plane, and the rotation axis of the crystal; 

Ó denotes the angle between the rotation axis of the crystal block and that of the 
crystal (the y-axis). 

The distribution function of the blocks, as computed from eqs. (14), (15), (16), 
(17), (18) or (22), (23), (24), results from the superposition of dispersion functions 
of blocks due to the different, variously oriented systems of dislocation lines. The 
distribution function of blocks due to a system of parallel edge dislocation lines 
is shrinked in the ratio cos ®, The angle ®, is the one subtended by the direction 


Wr 
A^ 


Fig. 8. 


of the dislocation lines and the axis of rotation of the crystal. The distribution func- - 
tion of those due to screw dislocations is shrinked in the ratio cos ®,. The angle Ø, 
is the one subtended by the normal to the plane of sliding containing the screw dis- 
location lines, and the rotation axis of the crystal. 


9. Computation of Dislocation Density for Various Forms of the Block 
Distribution Function 


If the block distribution function resulting from computation is of the Gauss 
curve form, then the mean angle between two neighbouring blocks can be computed 
from the formula (Gay, Hirsh Kelly 1953) $ i 


109 +060 
f Í eb x eby Ix — y| dx dy FR i 
on ack ooh eal re CE (2) 57098 x 


Í f et*&’+y?) dx dy 
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wherein fi is the half-width of the Gauss curve. Hence, by eq.(1), the density of dis- 
locations at the block boundary amounts to 


0.479 
Übound — bt P (35) 


By eq. (2), the total dislocation density amounts to 


. [0.479 B\? 
Q =a (36) 


For a uniformly bent crystal, the block distribution curve f(x) has the form a rectangle. 
The mean angle between neighbouring blocks is then 
2at 

Tig SD 
with 2a denoting the width of the function f(x), 
t — the block length, and 
L — the length of the crystal surface scanned by the X-ray beam. 
Assuming one dislocation to occur per inter-block boundary, the mean angle between 
the blocks is 

.. 2ah 


wherein À is the distance between the dislocations. Hence, the dislocation density 


in a uniformly bent crystal (Gay. Hirsh, Kelly 1953) amounts to 
2a 
e 7 (39) 


In practice, intermediate cases may occure, when, in addition to bending, there 
is dispersion of blocks in the crystal. It seems reasonable that, now, the resulting 
block distribution function f(x) should be described by the folding of two functions, 
namely, the Gauss function fj(x) and a rectangular shape function /,(x). If 


fila) = er” 
Jew) zc for z= as Pla) = 0 for el c>0 


then f(x) assumes the form 


+a a—x atx 
fi) = [core dy m efiet dete f.e" de (40) 
8 ò T] j 


Hence, p and a can be computed. 
The mean angle between contiguous blocks is now 


2\% 1 Lat 
À SMSI S I mae 1 
; a (2) itt T (41) 
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for t=h 


SAU : 
a= (2) Me > ar (42) 


Hence, the density at the block boundary is 


2 LER 
uw en SEEN te: 43 
Übound (2) pbt as Lb (43) 
The total density of dislocations amounts to 
lx 20 A2 1 MAE 
0= at tb a Td (n) 


Should the crystal contain no more than one system of edge dislocation lines 
parallel to its roiation axis, with Burgers’ vector lying in the plane of reflexion, the 
density as computed from eqs. (35), (36), (39,) (43) and (44) would be numerically 
equal to the number of dislocation lines intersecting a unit area of the surface, the 
latter being perpendicular to the rotation axis of the crystal. 

Should the crystal contain but one system of screw dislocstion lines lying in the 
plane of sliding, the latter being perpendicular to the rotation axis of the crystal, with 
Burgers’ vector perpendicular to the reflecting surface, the density as computed 
form eqs. (35), (36), (39), (43) and (44) would be numerically equal to the number 
of dislocation lines intersecting the reflecting surface. 

Nye (1953) proved that, for a crystal uniforly bent by a system of edge dislocation 
lines parallel to the bent surface, the curvature of the latter is given by the formula 


o = nb cos € (45) 


n denoting the number of dislocation lines intersecting a unit surface perpendiculary, 
b — the module of Burgers’ vector, and 
& — the angle between the bent surface and Burgers’ vector. 

In reality, a crystal exhibits a great variety of differently orientated systems of 


dislocation lines. By eqs. (33) and (35), the density as computed from eq. (39) can 
be written in the form of sums: 


pc X N,, cos e, cos ®,, + > > N, cos Ep, cos BD, 
r=1iel 

N,,is the number of edge dislocation lines subtending the angle @,, with the rotation 
axis of the crystal and whose Burgers vector subtends the angle e, with the plane 
of reflexion. ; 

N, is the number of screw dislocation lines in the plane of sliding whose normal 
subtends the angle ®,,, with the rotation axis of the crystal, and presenting a Burgers 
vector that subtends the angle &,, with the normal to the plane of reflexion. 
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The density of eq (35) can be representend in the form of sums, as follows: 


ft ns s 
bound = 22 i Na cos D, T > D NA cos C cos Dair 
‚v=1 


y-11-1 
The density as computed from eq. (36) asusmes the form 


my 


e = [X (4 Na) cos 5, + S >> (N,,)” cos é,, cos Dy)? 


v=1 v=ll=ı 
The density of eq. (43) assumes the form 


Ne Ns 


Übound = >: N,, (4 (1 AR q,) = q, cos é,)) cos o, Hs Di > Nes cos Epp COS Dar 
v=1 


v=1 1-1 


wherein q, is a quantity depending on the excess of dislocations of one sing for edge 
dislocations presenting a Burgers’ vector that subtends the angle e, with the plane 
of reflexion. 


10. Conclusions 


The method proposed yields the distribution function of mosaic blocks in a crystal, 
thus yielding the dislocation density therein. 

1) In order to determine the block distribution function, it is necessary in exper- 
imental work that the dimensions of the X-ray tube focus, spectrograph slit and 
photometer slit should be chosen sufficiently small for the following conditions to 
be fulfilled: 


a) ap<a 
wherein & denotes the angle of dispersion of the mosaic blocks, and f — the angle 
of divergence of the beam incident on the crystal in the horizontal plane: 


b) tg p cos? ply <a 
with 9 denoting the angle subtended by the beam incident on the film at the point 
under consideration on the line, and the one perpendicular to the axis of rotation 
of the crystal (Fig. 5). Ag is the angle of divergence in the vertical plane within the 
height of the photometer slit. 

2) If experimental conditions were such that the immobile and the oscillating 
film were situated at one and the same distance from the crystal, then the block 
distribution function is yielded by computation from eqs. (14), (15), (16), (17) and (18). 
The formulas were derived on the assumption that the outer surface of the crystal 
alone reflects (absorption is very strong), and that inequalities of the surface within 
the limits set by the width of the beam can be neglected. 

3) If experimental conditions were such that the immobile film was situated 
at one half of the distance from the oscillating film to the crystal, the block distri- 
bution function is obtained by computation from eqs. (22), (23) and (24). These 
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formulas also hold if there is reflexion from interior crystal faces (in crystals of low 
absorption); they eliminate error due to line deformation form modus tne of the 
surface within the beam width. 

4) According to the form of the block distribution function obtained, the disloca- 
tion density at the boundary between the blocks is computed from eq. (36), (39) or (44). 

The experimental curves required for computing the distribution function of 
mosaic blocks were obtained by measurement on several germanium monocrystals by 
Dr J. Auleytner and the present author, the numerical computations being carried out 
at the Institute of Computing Devices of the Polish Academy of Sciences (Zaklad 
Aparatów matematycznych PAN). 

The author wishes to thank Professor Dr L. Sosnowski for his valuable hints 
and discussions, and Dr J. Auleytner for suggesting the subject, and for his helpful 
advice 


REFERENCES 


Auleytner, J., Acta phys. Polon., 17, 111 (1958). 
Auleytner, J., Acta phys. Polon., 18, 81 (1959a). 
Auleytner, J., Postępy Fizyki, 10, 77 (19596). 


Burgers, J. M., Proc. Phys. Soc. [London]. 52, 23 (1940). "RET 
Gay, P., Hirsh, P. B., Kelly, A., Acta Met., 1, 315 (1953). PAN : ^ 3 
Nye, J. F., Acta Met., 1, 153 (1953). : 
P wit xi 
| 4 t ». $ ox ze 
, i m } y zt nm Ri ar hy ; eat Reiss 
M p d x "uet ga 2 


"E LAE Lm $n Mg ERES BB LI An t 


Vol. XIX (1960) ACTA PHYSICA POLONICA Fasc. 4 


ON THE DEPENDENCE OF THE PHOTON-ELECTRON 
RATIO ON THE DISTANCE FROM THE AXIS OF EXTENSIVE AIR 
SHOWERS OF COSMIC RADIATION 


By A. OLES 
General Physics Laboratory, Academy of Mining and Metallurgy, Cracow 
(Received November 11, 1959) 


The transition curves were determined for two distances of the measuring tray from 
the detector of extensive air showers. The ratio of the number of photons to the number 
of electrons was found to increase as the mean distance from the axis of the showers 
recorded was augmented (from 20 to 100 m). For the range of distances applied, a ~ 50 % 
variation of the ratio resulted. The sign of the effect observed is in agreement with that 
predicted by theory of electron-photon cascade. 


1. Introduction 


The problem of the variation of the ratio of the number of photons to that 
of electrons (f/e) in an electron-photon cascade as a function of the distance from 
its axis was elaborated theoretically by Eyges and Fernbach [3] and Borsellino ]2]. 
Their computations point to an increase in the excess of the photons over the electrons 
as the distance from the cascade axis increases. Notwithstanding the fact that -the 
mean angle of deviation for particles of a given energy is larger for electrons than 
for photons, the circumstance that the latter are able to travel a long way through 
the atmosphere without loss of energy leads (for cascades of the age parameter s >0.6) 
to an increase in the ratio f/e with increase of the distances from the axis of the cascade. 
Only young and narrow cascades with a small number of particles can exhibit the 
inverse effect. For these reasons, the foregoing conclusion of the electron-photon 
cascade theory should hold also with respect to extensive air showers of cosmic ra- 
diation. 

Hitherto, measurements of the fje ratio were carried out with the measuring 
tray located in the immediate neighbourhood of the detector of extensive showers. 
Milone [6] was the only author to investigate the dependence of the //e ratio on the 
distance from the axis of the showers recorded. Milone carried out measurements 
for two positions of the measuring tray, the one at the detector, and another one 
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50 m therefrom. The result obtained by him was in either case f/e ~ 1. With respect 
to the small difference of the mean distances of either position from the axis of the 
showers, Milone was unable to detect the effect under consideration. 


2. Arrangement 


Measurements were carried out at Cracow on the terraces of the principal bulding 
of the Academy of Mining and Metallurgy. Fig. 1 shows the set of the trays of the 
Geiger-Müller counters. The detector (T) of extensive air showers consisted of three 
trays A, B, C of an area of S — 0.45 m? each. A fourth, measuring tray (C) remained 


/ N 
/ rn a teas E E — Set ff 
(Irt bg RO 100 m Aud 


Fig. 1. Diagram of device. ABC — detector of extensive air showers. I, II — two positions of measuring 
tray. 


in coincidence with the detector, and was situated directly beneath an absorbent 
of variable thickness of 0 to 10 mm Pb. Measurements were carried out for two positions 
of the measuring tray: position I — within the triangle determined by the detector, 
and II — at a distance of 100 m therefrom. In the former case the area of tray G was 
% = 0.28 m?, and in the latter & = 0.91 m2. Aluminium counters, of a wall thickness 
of 1 mm and a threshold energy of — 1 MeV for electron recorded, were used. 


3. Experimental Results 


The numbers of coincidences registered are assembled in Tables I and II. 


Table I 


Results of measurements with tray G at position I. 


abite 
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Table Il 


Results of measurements with tray G at position II. 


Fourfold coincidences ABCG can by expessed by formula 
G() = k (le) (1 — oP RO*S) FD dy (1) 
0 


wherein ka~ + is the differential density spectrum of the showers recorded, S — 
the area of the detector tray, r = z %/S with z denoting the ratio of the mean densities 
of the particles over the measuring tray to that over the detector, and X — the area 
of the measuring tray G; R(t) is given by 


R (t) = P, (t) + fle P, (t) (2) 


wherein P, (t) and P,(t) are the probabilities that, an electron or photon respectively 
falling on the absorber of a given thickness ¢ will produce beneath it at least one 
secondary electron which operates the counter tray covered with that absorber; fje 
is the photon/electron ratio required, determined for energies of particles of either 
kind exceeding the threshold. energy of the electron recording device. 

By eq. (1), R(t) can be determined from the experimental data. The transition 
curves obtained are shown in Fig. 2. These served subsequently for finding the f/e 
ratio. In the computations (wherein eq. (2) was used), the experimental value of 
P, (t) was taken from the paper by Massalski and Oleś [5], and P, (t) was assumed 
as computed by Arley [1]. The following values resulted for the //e ratio: 


with the measuring tray at position I — (f/e); = 0.95+0.10 
with the measuring tray at position II — (f/e);y; = 1.404-0.15 


In an earlier paper [5], photons (especially of energy lower than the threshold energy 
of the electron recording device) were found to affect the coincidence frequency 
of the tray when not shielded by the absorbent. The effect involves an increase in 
the ABCD coincidence frequency for the uncovered tray, leading to a lowering of the 
maximum of the transition curve and yielding a lower value of the fje ratio. If the 
respective correction is accounted for, the values for (f/e); and (f]e)j; are raised by 
about 20%. — j 
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The mean distance of the measuring tray at position I from the axis of the showers 
recorded amounts to about 20 m. The analogous distance was evaluated also for po- 
sition II assuming the lateral distribution function for the particle density according 
to Nishimura and Kamata [4]. Computation yields for the measuring tray at position 
II a mean distance of ~ 100 m from the axis of the showers recorded. 


Ay I 


7,0 if 
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Fig. 2. Experimental transition curves. Curve I was obtained with tray situated within extensive air 
shower detector. Curve II — with tray 100 m distant from detector. 


4. Conclusions 


The experimental variation of the ratio of the numbers of photons to the electrons 
in extensive air showers of cosmic radiation as a function of the distance from the 
axis of the showers recorded exhibits a trend in agreement with predictions of the 
electron-photon cascade theory. 

i According to experimental data presented in Greisen’s article [4] the lateral 
distribution of electron density in extensive air showers at sea level is in agreement 
with the distribution for an electron-photon cascade of the age parameter s = 1.2. 
Similarly, the lateral distribution of photon density of extensive air showers be ex- 
pected to resemble the one in a single electron-photon cascade of the above value 
of the parameter s. Hence in extensive showers, the value of f/e should be expected 
to increase as the distance from the axis augments (as is the case in an electron-photon 
cascade of s — 1.2). > 

With the device described above, a ratio of f/e ~ 1 is obtained for the position 
of the measuring tray within extensive air shower detector. This, together with the 
increase of the fje ratio as the distance of the measuring tray from the detector aug- 
ments, leads to the conclusion that the photon/electron ratio throughout the entire 
shower exceeds 1. 
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J The author wish to thank Professor Dr M. Miesowicz and Docent Dr J. Massalski 
for their valuable advice and discussions; Eng. T. Owsiak and the technical team 
of our Institute for their help in the construction of the measuring device. 
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ON THE STATE OF A FERMI-SYSTEM WITH 
CORRELATION OF PAIRS OF PARTICLES WITH PARALLEL 
SPINS. I. GROUND — STATE. 


By ZYGMUNT GALASIEWICZI 
Joint Institute for Nuclear Research. — Dubna, USSR 
(Received December 9, 1959) 


The solution of a compensation equation of dangerous graphs, for a system of 
interacting Fermi-particles, as an odd function of wave vectors is investigated. It is shown 
that this solution leads to the anomalous state the energy of which is smaller than the 
energy of the normal state. The energy difference is equal to the binding energy of the 
pairs of fermions with parallel spin-moments. This anomalous state can exist in the pres- 
ence of attractive Coulomb interactions. The elementary excitations are separated from 
the ground state by a gap the width of which depends on the direction of the wave 
vector. For the direction perpendicular to the spin direction, energy gap is zero and 
therefore the state cannot be superconducting. The binding energy of the pairs and the 
width of the energy gap are of the same order as in the theory based on the creation of 
pairs with antiparallel spin-moments, where in the expansion of the interaction term 
in spherical harmonics, the first coefficient is more than three times greater than the 
zero coefficient. In the absence of the magnetic field the total spin of a system is zero. 
In the presence of the magnetic field the total spin differs from zero. The paramagnetic 
susceptibility is calculated. 


l]. Introduction 


Consider a dynamical system of Fermi-particles with the Hamiltonian 


BS DVL Gf) ap taper i Z Ui f ias fi) af af ap ap, (1) 
shes ps fif: 


TOS ys DOS f) ACE b 


where J is the Hamiltonian of a particle, U — the interaction energy, A — the chemical 
potential, ap, af the Fermi-amplitudes and f is a set of indices characterizing one- 
-particle states. In our case f = (p, o,), where p is a wave vector and a, a spin index. 


1) Permanent adress: Testitute of Theoretical Physics, Univesity of Wroclaw, Wrocław, Poland 
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Similarly as in [1] (and also in [2] and [3]), we transform the Hamiltonian (1) 
by means of the transformation 


dy = Dj (Uy, o copa). (2) 


To secure the canonical character of the transformation (2) the functions {u, v} 
must be connected by the orthonormality relations 


2j {us Up, ar Ufy DA Er vM = y 
x tus, Up, um Up, Up) = 0 (3) 


We find the functions {u, v) from the additional equations obtained from the 
compensation principle of dangerous graphs (see [2]) 


<a, a, H> 9 =0 (4) 


(the expectation value corresponds to the vaccum state in the &-represention) 
which are equivalent to the following equations for the functions F(f, f’), (f, f^): 


2i (8 6 f) 9 C05) + 808 PAY + SQ - 
= 2i QA) SG) + FA SUA = 9; (5) 
RS 
FAP) = 2; Op, ps P(A f) OFS) = us vp, (6) 
and | 
S (fifa) B U (fs fas fas fi) D (f,. fa) 
RIVERI) t GL SPS) “Uff pp Ege sf) c 


In order to represent the functions F, Ø in the form (6) these functions must 
satisfy the additional subsidiary conditions (see [2]) 


FEM = XU FUA 00.90.50) Ga) 
DEAN OLA + FAO LA}=0 (8 
From formula (1) of paper [2] we know that 
T(f.f) = (Ep) — ) (f — f) = elp) (f£ — f^. 
U (fi. fui fa fi) o Be Pisa EQ Obs toby Bir BO) d = a (oe o oe 


J (Py: Pai Pa» Pi) a J (Pay Pis Pr» Do) = J(— Py — Pa; u). (9) 
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The solution of equation (5) with the subsidiary conditions (8a, b) may be put 
in the form 


F(f. f^) = F(p) 6(f—f),® (ff) = 9(f) 8(f +f"), 
DE ot) (p, —) = o (p). (10) 


The functions F(p) and ®( p) are assumed to be real, invariant under the transform- 
ation of momentum reflection, and to satisfy the following equations 


D ANIM met 4; 
23 (p) ®(p) + ROE =D Pp) DP (p) 0 
F(p) = F?(p) + 6? (p). (11) 
where 
&(p) = e(p) + 7 [2/(p. p; pp) J(psp'sp,p)]F(p). . (12) 


The solution of this equations leads to the usual formulae of the theory of super- 
conductivity. 


2. The solution of the compensiation relation in the case ®(p, o,) is an odd 
function of p 


From the form equations (11) we see that a solution of this equations exist also 
in the case when (p, o,) is an odd function of p. Instead of (10), we put therefore 


F(f.f') = F(p)ó(f—f'). DS) = 9(f)8(p4-p)9(o— o). 
F(p) = F(—p). 9(p.o) = —9(p.—o). 6(p.o) =—®(—p,o,). (13) 
The functions of the form 


® (Dp, o.) = 9,0, DS.) PP) = 0, 0089 (D), p =(—B) (13) 


satisfy these conditions, where e = pjp. So, (Sz) is a spin function (S,, (o,) = 1, 
Sa, (s,) = 0 for o,#s,, 6, S = 0,S, 6,= +1), and ? is the angle between the wave 
vector and the quantization axis of the electron spins. Substituting (13) in (5) and 
(8a, b), we get for the functions F (p). D ( p. 9,) the equations (11) (from now on we 
shall denote the spin index simply by c). We define 


C. o) — — 3 MU ioe tat A: (us 


C(p, o) is an odd function of p and c 


C(p, o) = — G(— . o), C(p. o) = C(p, —o), C?(p, 0) = (217). (15) 
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From equations (11) we obtain 


(p A EATR Fp = 41 ae. 
POO oC) xq 7 20 VY Eat») 
hence 
1 A las 
reli nenn: (17) 
We take 
KB — Bs — BBY) = IP — Bl) = & 4G. P) P, (cos y) (18) 


where y is the angle between the vectors p and p' 
cos y = cos Ê cos # + cos g sin à son 0 (19) 


9, 9 are the angles between the vectors p, p’ and the polar axis, and g is the angle 
round the polar axis. According to the addition theorem, we have 


P, (cos y) = P, (cos 9) P, (cos 9) + 2 » e P5, (cos 0) Pp (cos 9’) cos mg. į 
mt , (20) 


Substituting (20) into (17) and writing the result in integral form, we obtain 
EO : 1 qme 


, n Cp, cos 0', a) p? dp’ ; 
= a De (cos 8) [J mown (cos 8") [Gres P ET sind’ dd’. 


par Fir 3 
d [ 


The solution of (21) may b put in the form Tu ohhh: "d 


19 


NS una ds E Des Fn 


= 
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where 


= dn 1 dn DE at 
un, 04 


(we have to replace the interaction of electrons with sound quanta by an effective 
electron-electron interaction, which is different from zero in a narrow layer near 
the Fermi surface Ep + c). In order that V — 0 for o, > 0 the effective interaction 
must be attractive (J;(pp) < 0). From the solution of (23), we get the function C(p, o) 
and consequently the functions F(p) and (p, o). Moreover, we can obtain another 
solution of (11) by putting (p, c) = 0 and F(p) = 1 for p < pp and F(p) = 0 for 
P > Pr: 

We examine now the special form assumed by the transformation (2) when 
the functions F and @ are given by (13). Taking into account (6) and (13’), we see 
that 


> 


uy, = u(B) 5(f—»), vp = v(— P, 0) ôP +) 6(0 — s), f= (P, 0), v = 3 


v(p, 0) = — v(—p. o), o(p, o) = — w(p, — o), wp. o) = a cos 9 w( p), w(p) = w(—p) 
(25) 

and 
F(p) = vp, o) = cos? 9 w? (p), ®(p, o) = u(p) v (p, o). (26) 


Therefore the transformation (2) has now the form 


az, = u(p) &zg— V (P, o) &*5,, a-zo = u (P) a3, + v(p,o)az,. (27) 


If by means of (27) we transform the Hamiltonian 


E OER 
H= M &(p) a$« Opa + 37 F. AE J (Bas Bas Bas Bi) Gyo, Gyo, Tro Azo, (28) 
bo Pob Pis 
(see (1) and (9)), we see that we must write the compensation equation for the dan- 


gerous diagrams for the pairs of fermions with parallel spins and antiparallel momenta. 
For the functions u(p), v(p, o) we get a first solution for the normal state, namely 


: Ri LEG) SE " £ 0, E(p) > E 
u(p) =O¢(P) = Is m 2 go P(B, o) — 8s. o) an Er cE, 09 


(where Op (p, o) = —9;, (—p, 0), ...), and a second solution by taking into account 
(26) and (11), namely 
> l ô (p) 2(3 1 t ô (p) | 
2(5) = — TER ARTEN E IT 
up 2 fı ii xl "Mis 2 V V? cos? 9 + 3* (p) 


é (30) 
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This solution is of the type obtained by Bogolubov [1] but it depends on the 
direction of the vector p and for 9 = 7/2 it goes over into the solution for the normal 


state. 


3. The energy difference for two states and. the energy gap 


Now we find the mean energy in the vacuum state in a-representation and 
evaluate the difference for solutions (29) and (30). The general formula for the mean 


value is 
E= <H >= WIEN 
IA 
dD Ulosa © Gale) OF pala) (31) 


ER 


Taking into account (9) and (13), we obtain 


<H>, = Y teo) all FD) +7 ȘI OE: PE) DB, o) Do) 
5 


pb 
& C? ( p) 
=2 AMD FG — 5) LIESS (32) 
p p 


The energy difference for the two solutions may be found with the help of the 
identity 
i.c) — OF (B)= 
29e], ——— ap | — Fe Ia iQ) oen 
VV? cos? 0 + 3?( p) 2 V V? cos? 0 -- 3?( p) 


Finally we get 


AE 2 à POOR 1 cQ 
— = — ) [v2 ( ,0)— 9 )I- = -A 
V pi VEM TW Legs 
= = FF dho? eye. (34) 


From this it is clear that the solution (30) leads to an energy state which lies lower 
than the normal state. This is also the anomalous state connected with the pairs of 
fermions with parallel spin-moments. 

Let us compare the ground state for the case of pairs of fermions with antiparallel 
spins with the ground state for the case of pairs with parallel spins. We see that the 
binding energy for the pairs of electrons with parallel spin-moments is proportional 
to [exp (—2/0,)]®, while the binding energy for the pairs of electrons with antiparallel 
Spin-moments is proportional to exp (—2/oj). When the case Ji(p, p) 2.3 Jo(p, p’) 
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is not realized the ground state of the system for pairs of fermions with antiparallel 
spin-moments has « lower energy than the state for pairs with parallel spin-moments. 

We shall now find the formula for the elementary excitation in the supercon- 
ducting state. The general formula is given by [3] 


<a, |H at >o = Y a f) (u up, — Uh Oy) + 
fF 
+4 23 Uh fos fif) IP (fa. /) (us v, — v ug) + 


fofas fifi 
4o (fa. Ar) (Up, Up» — Uy, y Up, ,)]- (35) 


By the way we remark that (35) is identical with the formula for the energy of 
elementary excitations given by Landau in [4]: 


00) == (36) 


Here E is the diagonal matrix element of the Hamiltonian (1) transformed by means 
of transformation (2) of [3] when the expectation value of æ} a, is not zero but n,. 
In this case the mean energy is given by a formula of the form (31) but the functions 
F and ® are now of the form (see [2]) 


D( fis J2) = > tur, vy (1 zn) Up Uy n,} , 
F(fi. fa) T >} {vrs Upe (1 = ny) I Uy y Ug v yj (37) 


For T = 0, that is n, = 0, formulae (37) fie over into (6). Taking also into account 
that in (31) F and p are functions of n,; we get 


5 = BED nip mn rn EEA 


Ony 2 
Lj à 
1 "vm ee / F1 1 Ot 6 (fi f2) 

ED URN en DFS + DA SEA | (ae 

fof» fifi E 
| Finding DNE from (37), we see that 

én,’ On, 
(3) = «a,|H|a; >, = Q (v). (39) 
ony ny=0 


Introducing (13) into (35), we obtain 


O(P) = (P) l^ B) (P. u 0) D Y. — 5: - P B») u (9). 
2 
P (40) 
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From (11) and (26), we get 


O (p) = VW? cos? 9 + 32(p) (41) 


On the Fermi-surface 
Q (p, cos 9) = |cos 9| 2 e** e^? = A (42) 


Thus the fermion excitations in the anomalous state are separated from the 
energy of the ground state by the gap A, depending on the direction of the vector D 
Since this energy gap can be equal to zero, no current-carrying state of the system 
will exist with respect to weak pertubations. This means that the anomalous state 
considered by us is not superconducting. 


4. The influence of the electron-electron Coulomb repulsion 


Let us now examine the combined effect of the electron-phonon interaction 
and the Coulomb interaction between the electrons (see [1] and [5]). Similarly as in 
the case of the electron-phonon interaction (which is essential only near the Fermi- 
-surface, Ep+w), we replace the Coulomb interaction by a model interaction. Because 
of the screening of the Coulomb interaction, we replace it by a constant attraction 
of electrons (essential only near the Fermi-surface F pw, c, > w). Thus we put 


A(t) = fas (3:8) + Ac( 3» lab]a| «e R 
J1 ` = = 
pie 6810.) 2 Se lal or |a’| >w n (43) 


We obtain now equations analogous to (23), from which we can determine PO 
in the two intervals: Meer w < |a| < o 


G uf 00) pg i 


MG puta ` EN. i e 


En 4 
Fo dM FRAIS W ors 


7 LS B PL 1 1 25 » d +i Le x1 
PEN Tome an Li X r * - Ee. "us 


S4 
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Hence 


1 20 l ) 

Zee 5 Pıc eu 

3 (m M EU | pean : a Or ; (47) 
Thus the condition of the appearing of a supereonducting state the for case of pairs 
of fermions with antiparallel spin-moments has an analogous form to that in [1]. 
It reads 


Bi ph s sg (48) 
l + X In 1 | 
3 01C PR 


5. Total spin of a system and paramagnetic susceptibility 


The spin moment for unit volume is given by the formula 
1 5 
S D = = 49 
y à, FG. - FG. (49) 
P 


We pass over to the operators œ and find (49) for the state C, (vacuum state in the « 
representation) 


-pY £F 0-FG.-). . 60) 
3 ; 


In the absence of a magnetic field F(p, +) = F(p, —) and therefore S* = 0. 

Consider now a system of Fermi-particles in the presence of an external, weak, 
constant magnetic field 6#, applied along the z-axis. In this case a new term appears 
in the Hamiltonian (28), namely 


ELPL LIE 6) 


where +e is the charge of the eas and m. its mass. 
The compensation equations lead to the formulae 


N e 
Aa ale ô (P) — — ot 
Hoy (P+) 5 5 ee ern; 
en 5 oe] eco 


y à(p) + — 6% 
hy (b. = 5 1 - er ee 1 
|o Foe] +e. 


' G(p, +) = C(p, —) ~ CD. (52) 
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Hence 


qx 
142 m 


sig : d} dx. (53) 
—o —1 | 


On perfoming this integration, we get 


d 
Soy = TE — ÒH [1 — exp (2/3 — 6/0,)] (54) 


For the paramagnetic susceptibility we obtain 


2e? dn 


t= a gp li 9 (2/3 — 6/o.)]. (55) 


From (54) we see that in the absence of a magnetic field the total spin of a system 
equals zero. In the presence of a magnetic field the electron-electron interaction 
cannot be taken into account by means of perturbation theory since we have ob- 
tained a non-analytical dependence of (54) and (55) for o, ~ 0. 

In this paper we have considered formally the solution of the compensation 
equations as an odd function of p. It is understood that from the point of view of 
transformation (2) this leads to a state of the system with pairs of fermions with 
parallel spins. In order to interpret the physical picture of this state one has to examine 
the thermodynamics of this anomalous state. 

It is a pleasure for the author to express his gratitude to Professor N. N. Bogo- 
lubov for suggesting the problem and helpful advice. The author also wishes to thank 
Dr D. V. Shirkov and Dr V. V. Tolmachev for valuable discussions. 
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AN ELEMENTARY DERIVATION OF THE FORMULAE FOR 
MULTIPOLE RADIATION 


By BRONISŁAW SREDNIAWA 


Institute of Theoretical Physics, Jagellonian University, Cracow 
(Received December 31, 1959) 


Expressions for the field intensities of the different electromagnetic multipole radia- 
tions produced by electric charges and currents as well as by electric and magnetic pola- 
rizations varying harmonically in time are deduced by developing the retarded potentials 
into series. The formulae obtained are direct generalizations of the Rubinowicz-Som- 
merfeld formulae for dipole and quadrupole radiations and are valid in the wave zone. 
The general formulae for the field intensities are also given in tensor form. 


Introduction 


In recent years the theory of multipole radiation has been extensively developed, 
particularly for the needs of nuclear physics. The results of these investigations 
are presented e. g. in the books of Blatt and Weisskopf (1952, chap. 8), Rose (1955) 
and Hamilton (1959, chap. 1), where also the literature of the subject is given. Since 
the methods of modern theory of multipole radiation are rather complicated I have 
tried to get the formulae for the field intensities of the electromagnetic field in àn 
elementary and straightforward way by simplifying and generalizing a method 
of Rubinowicz and Sommerfeld (Rubinowicz and Blaton 1932, Sommerfeld 1939, 
chap. 1). We assume that the electromagnetic radiation is produced by electric charges 
and currents as well as by electric and magnetic polarizations varying harmonically 
in time. 

We begin with the decomposition of the electromagnetic field into a „charge- 
-dependent" part and a „magnetization-dependent” part. For the first part of the 
field the retarded potentials are developed into series and decomposed again into 
two parts representing separately electric and magnetic multipole radiations. After 
the expressions for the field intensities of the „charge-dependent” part have been 
found, the field intensities for the „magnetization-dependent” part may be found 
without further calculations. No explicit use of the theory of groups and spherical 
harmonics is made. The formulae obtained are valid in the wave zone and are direct 
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generalizations of the Rubinowiez-Sommerfeld formulae for dipole and quadrupole 
radiations. Formulae for the electric and magnetic intensities of multipole radiation 
are also given in tensor form and are seen to be direct generalizations of the 
formulae for the lowest multipole radiations. 


§ 1. Decomposition of the electromagnetic field into „charge-dependent” 
and ,magnetization-dependent" parts 


Let the charge density o, current density jj, electric and magnetic polarization 


densities P and M be given in a spatial region V, outside of which all these quantities 
vanish. 
The Maxwell equations (in rationalized units) are then 


rot E+ H=—— M, div H = — div M 
les |; rae aide da 7S £ her 
rot H — E x (jp +P), div E = 0, — div P. (1) 
From these equations follows the law of conservation of electric charge 


9, o; + div jp = 0. 


We decompose each of the fields E and H into two parts: the , charge-dependent” 
part E,, H, and the , magnetization-dependent" part Ep, H,: 


satisfying the following equations 


IS: 


12% + : a 
rot H, — a E, — j, div E, = 0, "m (3) 
where 
ji P=j, o—divP =o, (4) 
and 


as jhe en 


1 
rot H, — > E, = 0, div E, = 0. (5) 
From (4) we see that the following relation is satisfied 


9,0 + divj — 0. (6) 
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§ 2. Retarded potentials and the electromagnetic field for the charge-dependent part 


The vectors E, and H, may be derived from the retarded potentials 


se. eg. f Bar, a(P, i) = z= 13 gy. (1) 


c Tpo 
where rpo is the distance between P and the source point Q and the expressions in 


brackets have to be taken at the time t—rpo/c. We assume that o, ... M oscillate harmo- 
nically in time: 


== ple 4-70" ef) Zee. (8) 

where g etc. are time independent. Introducing k = 27/4 = w/c we may write 
o= get + orc (8') 
Now we confine ourselves to the wave zone. We choose a point O so near to Q that 


> > 
r< R where r = OQ, R = OP (and R/R = n). It is well known (see e. g. Sommer- 
feld 1939, Vol. II, p. 61, Formula (12)) that the retarded potentials take then the 


form 
9 eik(R—ct) I ( Q en aS 9 


g-ikr:n — 3 E e (r . n)! 


Developing 


and denoting 


ó- 5 eik(R-ct) (10) 

we get 
yo 0X omn dr tere (11) 
a- S Yim nay +ere (12) 


Writing equation (6) in the form 
div j = ick o (13) 
we see at once that in our case the total charge vanishes: 


f odV =0 (14) 
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We shall decompose now (in a manner similar to that of Fierz 1949) each term 
of the series for a into a part which may be expressed with the help of o and another 
part independent of g. To this purpose we start from the obvious relation 


[ Bav=-— [rdvBav (15) 
valid for each B vanishing outside the finite region V. Substituting in (15) 
B=(n-r)'j, 1=1,2,... (16) 
and taking into account (13), we get from (15) and (16) 
s — ick | (nn! or dv = f (n:r) (m v) j + Un: j) r} dV (17) 


Now we may write the /-th term of a in the form 


= at fe ry dy =- CP Ics f e rn j) j + (n: j) v] aV + 


cfe ry} m 7)j — (nj) navy. (18) 


Taking into account (17) and putting n:r = r 
for d: 


= (— iyi 1 ch ORT, 
a=® [era toy EM [ exdxnav rc 
1=0 1=0 ke 


(19) 


„» we get finally the following series 


8 3. The charge-dependent field of the mu‘sipole radiation 
a. Electric 2'-pole radiation > 


Let us take together the successive terms for y and for a from (11) and the first 
part of (19): 


mI. 
wo E faar teo, PATON (20) 


eT - 
ap = 9S [tend tere (21) 


These potentials satisfy the Lorentz condition in the wave zone. Thank to the relations - 


e ex 
E, = — grad p — — å, H, = rota (22) 


ngo o 
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we get from (20) and (21) the field intensities 


Eo! — (Gana N 1-1 i 
(i — — TI ri-* om dV+e-e. (23) 
Hl pp exta, - e)l 

( ð " rn-InxotdV 4- c.c. — nx ES (24) 


where r, is the component of r perpendicular to m. 
These fields satisfy the symmetry relations 


EQ" (—R) = (—1'? EQ, H?(-R)-(-1'H9q —— Qs) 


characterizing the electric multipole radiation. 
Taking into account (4), we can finally write (23) and (24) in the form 


KS pct et (gO py dv (26) 
OL = age 5? 1 (or — div P) r, dV + c.c. 


Ho ep (27) 
(o — d dem Iin of — div P) dV + c.c. 


b. Magnetic 2'-pole radiation 


The remaining part of @ in (19) yields 


gy? ES 


jl (— ik) l 3 T] ` 
(m) l .C. zb 2 sss 28 
qua ar cm fs (r xj) x ndV + cc, l=1 (28) 


the Lorentz condition for 9” and al” is again satisfied and simple calculations lead 
from (3) to the following expression for the field intensities 


Emi= p — (kd ri nx {rx (jp — ick P)) dv CHE (29) 
KED 1)! n Jf i ze 
Hm ——o (= iy ttr ipti r-1 (rx( jj —- ickP)) dV + c.c. = nx EQ (30) 
cm (l Die " Jf GL: ] 


which satisfy the symmetry relations 
EP (—R) = (VER, HP (—R) = (VER) 6) 


characterising the magnetic multipole radiation. 
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c. Analogies 


Comparing formuale (23) and (24) on the one hand and (29) and (30) on the 
other hand, we see that the replacing in (23) and (24) 


of El, HE EN, 


l : 
d ey i A 
by Hi. E, (i+ Bas 
gives exactly formulae (27) and (28) for all multipole radiations. j 
84. The magnetization-dependent part of the multipole radiation | 


Upon replacing 
| E, H, M, — div M. 
| by (32) 
: —H, E, —j, e 


equations (5) change into (3) and we can write directly the expressions for the magnet- 
ization-dependent part of the multipole radiation without further calculations. 


| i a. Electric multipole radiation 

i |. From (29) and (30) we get | : roca E E 
nah lee. T 

E¢ ® rar jor (rx Ma sc ISIN 


Uie to GH i wo La 
RU +D fe ler ar a) 
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c. Total field 
We obtain the total electric or magnetic 2"pole radiation by superposition 
Eo! = Ee a EQ, Ho — He Je He 
Eo — Em! JJ Et, Ho em Ho! n. Hy! (37) 


From these formulae the energy and the Pointing vector for different multipole 
radiations may be easily calculated. 


$ 9. Multipole moments. Tensor form of the formulae for multipole radiation 


We define the time-independent multipole moments as follows: 


a. Electric multipole moments 
ee | Aj... Xi (0g — 9, P) dV 


2 mE - = 
id f Kia Th, Ez M T$ He jak (38) 


b. Magnetic multipole moments 


"HE IE en LU i Ural ick. F; x) m Gi, — d: P ak 


31-441 


a= | 6,9; M; dV (39) 


(All indices run form 1 to 3). These multipole moments are tensors with respect 
to linear spatial transformations. The multipole moments depending on time are then 


1o, gd --c*c. ete. (40) 


ie) ije Ay 


Introducing the completely antisymmetrical tensor &;,, taking account of the space- 
-time dependence of ®, and putting 


nj vee ng 
Ai, = a RIT’ 
I--1 
Qe» ES d N! 
Qr.. dtiti ego 


we see that we can write the formulae for the electric and magnetic multipole radiation 


in tensor form as follows: 
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HOY lái t eee 
a. Electric 2'-pole radiation (from is ), (33) and (27), (34)): 


7104-1) QS 
E09! = Aj... a? Qro Dow wes (l = as Qiu ARE OR 


c 


N2U+2) 
HO! = A, won [ote ik iki, — TIC Qu. J _R 


b. Magnetic 2! — pole radiation (from (29), (35) and (30), (36)) 


x 204-1) FU) 
Ll Ai. 3 Uu dk fihi, — Te: Ji,. a if, -2 1 
De vr A41) Sagrada) ss [s 1 
Ho = AL lr: Je tia m TE t), R | 
I 


These formulae are the generalizations of the well known formulae for the lowest 
multipole radiations. 


§ 6. Final remarks ^s 


a. We may get rid of div P and div M in formulae (26) and (27) by ang 
| > of the identity 


: } he MUT ` | a 
for a — r-!r, ar Uns ire and be -P or b=M. 


We ge then after S calculations instead. of ee and md `~ i d : 
‘Sat ; ! X 4s ENS 
penes 
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and similarly for a(?, af” differentiating the “Maxwell R” (in the exponent) instead 
of the “Coulomb R” (in the denominator). 


I wish to express my thanks to Professor Jan Weyssenhoff and to Miss A. Ko- 
walska for interesting discussions and remarks. 
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EINFLUSS DER KRISTALLGITTERBEGRENZUNGEN AUF 
DIE SPINWELLENRESONANZ IN EINEM FERRIMAGNETIKUM 


Von A. R. FERCHMIN 


Institut für Theoretische Physik der Adam Mickiewicz Uniwersität, Poznan 


( Eingegangen am 13 Januar 1960) 


Es wird die Spinwellenresonanz untersucht, welche von dem homogenen magnetischen 
Wechselfeld in einer dünnen ferrimagnetischen Schicht erzeugt wird. Man erhält eine 
Lösung von N halbklassischen Spinbewegungsgleichungen. Die Form dieser Lösung 
hängt davon ab, ob an beiden Enden einer linearen ferrimagnetischen Kette die Spin- 
momenten gleich oder ungleich sind. Es wird gezeigt, dass man ausser den gewöhnlichen 
Spinwellen mit einer periodischen Amplitude auch die aperiodische Spinwellen erregen 
kann. 


J. Einleitung 


Die den Kristall begrenzenden Flächen verursachen Spinwellenreflexionen und 
unter bestimmten Bedingungen können stehende Spinwellen entstehen. Die Spins 
an der Oberfläche befinden sich unter der Wirkung anderer Kräfte als die im Innern 
des Kristalls, und zwar wegen der fehlenden Nachbarn. Dieser Gedanke wurde zuerst 
von Kittel (1958) ausgesprochen, und er hat ihn zum Voraussagen der sogenannten 
Spinwellenresonanz im homogenen Hf Feld geführt. Zwischen dem halbklassischen Bild 
der Spinwellen und den mechanischen Wellen in Kristallen besteht eine Analogie. Das 
Problem der Wellenausbreitung in endlichen Kristallgittern ist schon lange bekannt 
(siehe z.B. Born and Huang 1954). Man sieht in endlichen Gittern von der sogenanntem 
Periodizitätsbedingung ab und setzt andere Randbedingungen voraus. Der Fall einer 
endlichen linearen Kette wurde auch von Wallis (1957) diskutiert, der u.a. einige 
Sonderlösungen, die er Oberflächenschwingungen nannte, gefunden hat. Kittel (1958) 
bediente sich der Bewegungsgleichung für den Spin in einer Form, die für stetige 
Medien gilt, d.h. wenn die Gitterkonstante a nach Null oder die Wellenlänge 
nach Unendlichkeit strebt. Die Behandlung eines Kristalls wie ein stetiges 
Medium, was bei grossen Kristallen begründet ist, kann in dünnen Schichten 
manche Diskontinuitäts- und Oberflicheneffekte (z. B. ungleiche Differenzen der 
sukzessiven Wellenzahlen der stehenden Wellen) verschleiern. Kittel nimmt ein 
Ferromagnetikum in Betracht; Orbach und Pincus (1959) und unabhängig von 
ihnen Cofta (1960a) haben den Fall eines antiferromagnetischen Stoffes behandelt. 
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Der allgemeine Fall eines ferrimagnetischen Stofles, in welchem die anderen Fälle 
mitenthalten sind, wurde von Cofta (1960b) diskutiert. Cofta nahm die von Kittel 
(1958) formulierten Randbedingungen an, und löste dann nur die Bewegungsleichun- 
gen für die inneren Spins. Auf diese Weise nahm er nur indirekt den Einfluss der 
Kristallenden in Betracht. Eine solche Betrachtungsweise ist aber vom mathemati- 
schen Standpunkt nicht begründet, da die Spinbewegungsgleichungen, wie wir weiter 
sehen werden, miteinander gekoppelt sind, so dass sie ein System von Gleichungen 
bilden. (Siehe auch Fox 1957 S. 5) Dieses System kann nur dann in einzelne Gleichun- 
gen aufgespaltet werden, wenn alle Gleichungen (bzw. Gleichungspaare) dieselbe 
Form haben, d. h. wenn man von den Randgleichungen absehen kann. Für sehr 
dünne Schichten kann aber eine solche Annahme eine ungenügende Annäherung 
sein. 

Wir wollen hier die oben erwähnten Annäherungen vermeiden und eine ferri- 
magnetische endliche lineare Kette untersuchen. Es wird auch der spezielle Fall einer 
ferromagnetischen Kette untersucht. 


II. Spinwellen in einer endlichen Kette 


A. Ungleiche Spins an den Enden 


Betrachten wir eine lineare Kette in der Richtung der z-Achse von äquidistant 
fixierten Spins von zwei Arten: S, und Sz, die alternierend gesetzt sind, so dass ein 
Atom mit dem Spinmoment S, zwei S, als nächste Nachbarn hat und umgekehrt. 
Die Bewegung eines Spins wird durch die bekannte halbklassische Gleichung 


aS 2J x 
eu Ñ Surety Six H (1) 


beschrieben, wo & die Summe der Nachbarspins und J den Austauschintegral zwischen 
Nachbarspins bedeutet. Die Spins an den Enden der Kette kónnen ungleich oder gleich 
sein. Wenn sowohl die Zahl der S, Spins, als auch die Zahl der S, Spins N beträgt, 
so ist der erste Fall realisiert. Die Spins befinden sich in einem zur z-Achse parallelen, 
statischen magnetischen Feld H,. Dieses Feld enthält auch das effektive Anisotropie- 
feld, das auf die inneren Spins wirkt und darum unterscheiden wir H,, und Hy. 
Wir folgen Kittel (1958) und nehmen an, dass in der z-Richtung das Oberflüchenani- 
sotropiefeld H,, auf den Spin S,,, und H,, auf den Sj; wirkt. Wenn man annimmt, 
dass S? = S, = const. und 5j = S, = const., findet man mit Hilfe der Formel (1) 
für die y-Komponenten der Spins: E 


9S vA petes x x 
am ee ELI: (Sar Sp — Spo Se) — YSaı (Hao + His) 
9S; 2;—1 2] x x 3 
ES = — | [2 Sazi- Sp — (Sp 2j-2 + Sp2j) Sel — y Sa aj-1 Heo 
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9S3»; DIC { 5 
3t Tags à [2 Sf 2j Sa — (Se2j-1 + Saaj41) Sp] — Y Sia Hgo 
OSB on PA PUT US x 
HT OPN Sa — Se2n-15p) — y Span (Hgo + Hap) (2) 


und ähnliche Systeme von 2N Gleichungen für S* und $*. 
Setzen wir 


x MEN iet 


, 


x = iot 
SB 2; = Ag, e 


Saj- = Pda, ee”, 
| Sa Awe: (3) 
so kónnen wir schreiben: 
(u + a) A, + 45 — 0 
As; + u Ag; a + Ag = 0 
Ag;_1 + V Ag; + Aag = 0 
An-ı + (w + b) Any = 0, (4) 
(2N Gleichungen) 
wo 
u = (o — 20,5 + Wq0)/Deq> V = (w — 20,4 — c)/0,5; 
a = (Wep — Vaa) Oea» b = (Wea — O,5)/ 0,5; 
Og = (J/h) Sas Og = (JIA) Sp, 
Wao = y Hag. Opon 4 Ho» Oe = y (Hot Ha) Qg — y (Apo ru His); 
Waa = Og — Vays Vag = (0g — Wp: (5) 
Das System von Gleichungen fiir die S*— Komponenten hat dieselbe Form (4) und 


das für die 5 Komponenten geltende System ist identisch erfüllt von den Lösungen (3) 
mit beliebigen Amplituden A. 


Die säkulare Gleichung für das Gleichungssystem (4) ist: 


¥ 


Dy (a, 6) = 0, (6) 
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Dina) = I E (7) 
Lees eel 
1 v+5|2N 


(2N markt an die Ordnung der Determinante). 
Wir folgen jetzt Wallis (1957) und Rutherford (1947; 1951) und führen einen Para- 


meter © ein, indem wir 


(uv) = — 2 cos O (8) 


seizen. Mit Hilfe der Rechenmethoden von Rutherford kann die Determinante (7) 


nun in der Form 


uv + bu + av 
— (uv)* sin © 


sin (2N — 1) O 


Dan (a, b) = sin o (9) 


sin 2N O + (ab — 1) 


ausgedrückt werden. 

Betrachten wir zuerst einen Sonderfall, wo die Oberflüchenanisotropie klein 
ist, d.h. œ, %w,. Gerade in dieser Annäherung hat Kittel (1958) den ferromagneti- 
schen Fall untersucht. Es folgt dann (s. Ausdruck 5) a = 6-1, ab = 1 und die säkulare 
Gleichung (6) wird in die Form 


sin 2N O Á 


(uv + bu + av) Soy Ve 0 ` (10) 


übergeführt. 
Die Lósungen dieser vereinfachten Gleichung findet man aus der Beziehung: 


sin2NO0 _ o 
sin 20 TH 
oder 
E TN EUR pea 
Die Lósungen von (11) sind: 
O,-—n(nN) lxmnsN-—I (13) 


wo n eine ganze Zahl ist. Aus der Formel (8) folgt mit (5) (ohne Annäherung w, < w,) 


a =o [14 Oe E (on tn) 
OE 


(cos? O + e? sin? ey | : (14) 
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wo 


a [Wee NN i (Mao cs 059)]? 


E 
[Wea + Dep — $ (Wao — @go)]? + 2weg (Wao — cgo) ' 


CO T Op a Ong 21360 59)- 


Die Vorzeichen + gelten fiir zwei Freqenzzweige, welche wir wegen der Analogie 
zu den elastischen Wellen in zweikomponentigen Kristallen, akustischen und optischen 
Zweig nennen werden. Diese Formel (14), die sich auch in der Weise 


Q0 -—3- ilo. Dep i (a0 go) |” T Aea Dep cos? 9j : (15) 


schreiben lasst, stimmt mit der Formel (9) von Cofta (1960b) überein, wenn wir den 
Parameter © durch k, a (k,— Wellenzahl) ersetzen. Der Ausdruck (14) gibt uns 
also eine physikalische Deutung des Parameters ©: er ist dem Spinwellenimpuls 
proportional. Wir werden ihn ferner kurz „Impuls“ nennen. Der Ausdruck (14) 
stellt also eine Dispersionsformel dar. Wenn wir in (15) S, = Sg bzw. S, = —S, 
setzen, so bekommen wir Formeln, die entsprechend für Ferromagnetika oder Anti- 
ferromagnetika gültig sind. Eine ausführliche Analyse der Formel (15) befindet 


sich in der Arbeit von Cofta (1960b). 


Als mögliche Lösungen von (12) findet man leicht: 
WD = Wy, (16) 
w = Oy = Wea + 0,5 + Op- (17) 


(In der benutzten Näherung können wir «X, = @g = Wo setzen). Wenn wir in (14) 
die N—1 Werte der O von (13) einsetzen, bekommen wir 2N—2 Werte der w, welche 
mit zwei Lósungen (16) und (17) die 2N Werte der Frequenz bilden. Die Lésung (17) 
liegt in dem ,verbotenen* Bereich zwischen dem akustischen und optischen Zweig 
der Frequenz. Der dieser Lósung entsprechende Impuls ist komplex und wird durch 
die Formel 


O = ln + iar sin hyp (q/2) (18) 
ausgedrückt, wo | 


q = (5,— Sp) (Sz Sp)? 


Wenn wir die Oberflächenanisotropie nicht vernachlässigen können, lässt sich 
die sákulare Gleichung nur graphisch lösen. In diesem Falle bekommen wir 2N 
Werte für ©. Sie sind nicht ganz genau gleich voneinander entfernt, wie «diejenigen 
im Ausdruck (13). Falls die Oberflächenanisotropie, wenn auch nicht vernachlässigbar, 
aber klein genug ist, bekommt man u.a. auch eine komplexe Lósung. 
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B. Gleiche Spins an den Enden 


Untersuchen wir jetzt eine lineare Kette mit gleichen S, — Spins an ihren Enden. 
N sei eine gerade Zahl. Die Kette besteht aus N+1 S,—Spins und N S,—Spins. 
Die Bewegungsgleichungen für die inneren Spins sind in (2) gegeben. Die Bewegungs- 
gleichungen für die Spins an den Enden haben aber die Form: 
955 2] 
en 


Sai Sp — Sa Sa) — y Set (Hao + Haa), 


9 Se aN41 


Dall x Pal x X 
orn E EN (Sa 2N+1 Sp — Span Se) — Y Sa anti (Hao + Has). (19) 


Setzen wir in diese Differentialgleichungen die Lósungen (3) ein. Unter Verwendung 
von (5) kónnen wir schreiben: 


(u + a) 4, + 4, —0 


Aj. X UA, ; 4 it A; =0 


RN VETERAN VS 


Agy + (u + a) Ao 1 —0 (20) 
(2N+1 Gleichungen) 
Wir führen nun, Wallis (1957) folgend, die symmetrischen Präzessionsradien 

Ag m [45.1 SE Ayn-j+9-1] 

Ay $ [45 + Ast Lj) (21a) 

mit 1 <j < (N/2), 
Ag = Aja (21b) 
mit j = (N2) + 1 


und die antisymmetrischen Präzessionsradien 


Aaj de [45-1 5: Ag v L2) 3: 
Ag; = 5 [43; — Ayn_j+y] (22) 
mit 1<j<(N/2) ein. Die .4* erfüllen ein System von N-+-1 linearen homogenen algebrai- 


schen Gleichungen. Die A erfüllen ein anderes System von N homogenen Gleichun- — 


_ gen. Beide sákulare Gleichungen für die symmetrischen und antisymetrischen Lósun- 
gen A* und A lassen sich in der Form: 


(u/v)* cos (N + 1) O — a cos NO = 0 (23) 
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(für symmetrische Prüzessionsradien) bzw. 


sin (N4-1) © , sin NO 
sin © a eru) snQ 0 (24) 


(für antisymmetrische Präzessionsradien) schreiben. Aus (23) und (24) erhält man 
mittels graphischen Methoden 2N--1 Werte von ©, und zwar N--1 Werte von © im 
Bereich 0 < © <z/2 (akustischer Zweig), und N Werte von © im Bereich (r/2) <O «cz 
(optischer Zweig). Man bediente sich hier ebenso wie Wallis (1957) des sogenannten 
»Schemas ausgestreckter Zonen“. Für kleine © aus dem akustischen Zweig gilt an- 
nähernd die Formel (13). Diese Annäherung ist um so besser, je kleiner O, und je 
langer die Kette ist. 
Für die symmetrische Wellen erhält man: 


cos (N — 2j -- 1) O 
cos NO i 


cos (N — 2j +2) 0 
cos NO - 


A = O(S 


Arere C(— 1/7! uoi (25) 
Es wird später bewiesen, dass die antisymmetrische Wellen durch das äussere Wechsel- 
feld nicht erregt werden. Deshalb brauchen wir sie nicht berücksichtigen. Die Form 
dieser Lösungen zeigt uns noch einmal, dass @ der Wellenzahl proportional ist. 

Falls S,, Sg > 0, d.h. die Spinmomente gleichgerichtet sind, und die.Kette ge- 
niigend lang ist, haben wir nur N reelle Wertee von O im akustischen Zweig und N—1 
reelle Werte von © im optischen Zweig. Die zwei fehlenden O-Werte sind komplex. 
Durch Substitution von © = (z/2)--iy in die Gleichung (23) für die symmetrischen 
Wellen, und © = (z/2)— iy in die Gleichung (24) für die antisymmetrischen Wellen 
erhült man Gleichungen, welche man graphisch lósen kann und auf diese Weise 
findet: man jene zwei komplexen Werte von ©. Die entsprechenden Werte von w 
(sie sind immer reell) liegen zwischen dem akustischen und dem optischen Zweig. 
Wir erhalten dann gedämpfte Wellen mit einer Amplitude, welche von den Enden 
der Kette zu ihrer Mitte abfällt. Wir werden sie „Oberflächenspinwellen“ nennen. 


III. Sonderfall: eine ferromagnetische Kette 


In diesem Falle haben wir S, = 55 = S, und (ujv) = 1. Die Wellenamplituden (25) 
haben also eine vereinfachte Form. Wir wollen sie mit den Lósungen Kittels (1958) 
vergleichen. In (25) haben wir aber den Anfangspunkt der Koordinaien in der Mitte 
der Kette, während er sich bei Kittel auf der Oberfläche befindet. Darum müssen 
wir den Cosinus anders ausdrücken; nämlich: 


cos (N — 2j + 1) O/cos NO = tg NO sin (2; — 1) O + cos (2; — 1) O. (26) 


Wenn O, = n(r/2N), dann ist tg NO für ungerade n sehr gross. In diesem Fall also 
spielen nur die Sinuswellen mit ungeraden n eine Rolle. Das stimmt mit der Behauptung 
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Kittels überein. Aus einer Analyse der graphisch ermittelten Werte von © ergibt 
sich, dass die Annäherung (13) um so weniger genau ist, je grösser das n und je kleiner 
das N ist. Deshalb können in sehr dünnen Schichten bei Resonanzen höherer Ordnung 


auch Cosinuswellen und Sinuswellen von geradem n erscheinen!. 


IV. Erregung der Spinwellen 


Nimmt man an, dass ausser dem statischen magnetischen Feld H, auf die Spins 
noch ein äusseres, zirkular polarisiertes Feld in einer zur z-Achse senkrechten Richtung 
einwirkt, so betragen die Komponenten des totalen Feldwektors H, welche auf einen 


S,9pin wirken: 
H, = hy exp (iwt), H, = ih, exp (iot), H, = Hao (bzw. H, = Hao + Haa) 


Einen analogen Ausdruck erhalten wir fiir einen S,-Spin. Nehmen wir zuerst eine 
Kette von 2N Atomen in Betracht. Die Bewegungsgleichungen für die Spins sind in 
diesem Falle nicht homogen. Wenn wir als Lósung den Ausdruck (3) mit 


Aja => A571 "id B, 
und Asi As Bs (27) 


ansetzen, finden wir 


(u + a) A, 4- EE C, 
Aja tu Ayj-ı a A5; = 
Asi "m Aj; ar Aja 229 
Ze quc by Aon = Cs (28) 
(2N Gleichungen) 
Wir benutzten hier folgende Bezeichnungen: 


Q5 = phy, 


m h Cg Weg 


[ee ZF D [w t i wo T 2 (Wea an Wes)|, 


h wwe 
Bg == aD) P [o Wego — 2 (Oea = ®ep)|, 


: 1) Hert Dr. H. Cofta zeigte dem Autor eine während der Bearbeitung dieser Arbeit noch nicht - 
veröffentlichte Arbeit von Pincus (1960). Seine Analyse der Resultate von Kittel (1958) ist der hier 
angegebenen Analyse ühnlich. Pincus geht von der Symmetrie des Hamiltonians aus und bekommt 
einem Ausdruck, welcher der Formel (26) analog ist. , E 
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h On Dey 


Ca => 2J D {ae [o e230 a 2 (Wea "E Wes) | = (eg (€ go == O0) - 


h Wh Oek i s 
Cs = S SEDI. {@up |0 — wao — 2 (Mea + ep)| + Wee (920 — @po)}: 


D = o? — o |o, + cgo + 2 (cg + €,5)] + Og Mp0 + 20.4 Oro + 20,559. (29) 


Das Glied B, oder B, im Ausdruck (27) hängt nicht von j ab, daher stellt es eine 
unendlich lange Welle dar. Das äussere Wechselfeld tritt in (28) nur im Ausdruck 
C, und C, in der ersten und letzten Gleichung auf. Dieses Feld erzeugt also Wellen 
endlicher Länge.nur auf der Oberfläche. Ihre Amplituden A’ erhält man aus (28): 


; —l . , sin 2N O sin (2N — 1) O 
Ay = i | 
: Dax (a. b) [e LZ sin O ^ sin O | i c) 
; MET "o dn(2N—2/2)0 — ,snN—3j4- 109 
Te 06 : E 
23-1 ~~ Dax (a, b) | Š [em sin O h sin O i 


, | sin (2j — 1) 0 ‚sin (27 — 2) © 
ib | m O aaoo: E Oo | 


je =) sin (2N — 27+ 1) O sin (2N — 2j) ©. 
Aa; Don (a, b) | x | sin O rule sin O * 
2j 0 sin (2j — 1) O 
+ Cs [ u[v)* T a A acon en ev 


ESTO MEN du uw xm her ^ og ch MEN 5 wq 94.9 "XAR E "s e eA fe 


F E sin2N@  snQN—1)0 
Aon = Dos: (a, b) IC. d Cs joy = On =e Alt (30) 


Nehmen wir jetzt an, dass die Kette aus 2N--1 Atomen besteht. Durch Anwendung 
der Formel (3) und (27) bekommt man für die inneren Spins die, Gleichungen (28), 
sowie die Gleichungen 


(u + a) 4, + 4; = Ca» 
Ay + (u + a) Aog = Ca (31) 


für den ersten bzw. letzten Spinpräzessionsradius. Wenn wir jetzt die symmetrisierten 
A'* und antisymmetrisierten A’~ Amplituden, welche durch die Formeln (21) und (22) 
mit gestrichelten A’s definiert sind, anwenden, bekommen wir neue Gleichungen. 
Die Gleichungen für 4’ sind homogen, so dass äussere Wechselfeld in diesen nicht 
auftritt. Auf diese Weise stellen wir fest ‚dass das äussere Wechselfeld die antisymmetri- 
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sche Wellen nicht erregt. Die Amplituden symmetrischer Wellen kónnen wir mittels 


der Formel 


e XA 
s z — (2j — 2] 6. 
Ana a cos NO — (ujv)? cos (N + 1) O uns ) 
A5; = Ca cos [N — (2j — 1)] O (32) 


a (v/u)* cos NO — cos (N + 1) O 


ausdrücken. Der Resonanznenner in (30) und (32) ist gleich Null, wenn O Werte 
annimmt, welche die sikularen Gleichungen (6) bzw. (23) erfüllen. Wenn diese 
Werte komplex sind, dann gibt es eine Resonanz auch für die aperiodischen ,,Ober- 
flàchenspinwellen*. 

Wir versuchen jetzt die Tatsache aufzuklüren, dass wir andere Ausdrücke für die 
Spinwellenamplituden erhalten, je nachdem, ob sich an den Enden der Kette gleiche 
oder ungleiche Spins befinden (vgl. Formel 32 und 30). Jede stehende Spinwelle im 
Innern des Kristalls besteht aus zwei fortlaufenden Spinwellen, die durch das áussere 
magnetische Wechselfeld auf seinen Begrenzungsflächen erzeugt werden. Wenn sich 
an beiden Enden der Kette gleiche Spins befinden, entstehen zwei gleiche Spinwellen, 
welche eine stehende Spinwelle mit einer Amplitude (32) ergeben. Wenn es an den 
Enden der Kette ungleiche Spins gibt, entstehen zwei Spinwellen mit ungleichen 
Amplituden und geben eine Spinwelle mit der Amplitude (30). 

In dieser Arbeit haben wir nur die Austauschwechselwirkungen zwischen den 
nüchsten Nachbaratomen .berücksichtigt. Es ist bekannt, dass in den Ferriten und 
in den antiferromagnetischen Substanzen eine wichtige Rolle die Wechselwirkung 
zwischen den zweiten Nachbarn spielt. Es scheint, dass zur Lósung dieses Problems 
man die Rechenmethoden von Rutherford (1951) anwenden kónnte, unter Berück- 
sichtigung sowohl der Austauschwechselwirkung zwischen den erstnächsten, wie auch 
zwischen den zweitnüchsten Nachbarn. 

Die Rechenmethoden von Rutherford kónnte man. ebenfalls zum Problem der 
Spinwellenresonanz in flachen Netzen und in einfachen dreidimensionalen Gittern 
anwenden. Kittel (1958) zeigte, dass die in den dünnen ferromagnetischen Schichten 
parallel zur Oberflüche laufenden Spinwellen keinen wesentlichen Einfluss auf den 
Verlauf der hier beschriebenen Erscheinung ausüben. Man kann erwarten, dass dies 
ebenfalls bei den Ferrimagnetika und Antiferromagnetika der Fall sein wird. Das 
wird wahrscheinlich der Grund sein, weshalb Cofta (1960 a, b), Orbach und Pincus 
(1959) sich auf die lineare Kette beschränken konnten. 

Herrn Prof. Dr. S. Szczeniowski danke ich für die Anregung zu dieser Arbeit, 
seine kritischen Hinweise und für das gefällige Durchlesen des 'Manuskriptes. Herrn 
Dr. H. Cofta danke ich für die Einführung in das Problem der halbklassischen Spin- 
wellentheorie und seine noch nicht veröffentlichten Arbeiten, die er mir gern zur 


Verfügung stellte, Herrn Dr. L. Kowalewski und Herrn L. Nowicki für nützliche 
Diskussion. | ] 
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An equation, invariant with respect to the 12-parameter CXC’ group, for a spin 
1/2 field interacting with an electromagnetic field is investigated. The groups € and C’ 
are treated symmetrically. Correspondence with Dirac’s equation and bilocal theory 
is disscused. We obtain the result that mass is an isovector, in a similar way as in Rayski’s 
generalization of Pais’ theory, and charge is an antisymmetrical isotensor of second 
rank. The charge and mass conjugation and the existence of a charge gauge group are 
investigated. 


Introduction 
In Dirac’s equation 
Dy TTE (Yu 9, yt m) s 0 (0.1) 


the mass m is treated as some untransformable quantity. Thus we conclude that 
the groups of continuous transformations, permitted by (0.1), are 


C: y = eut ay 8, = a,, 2, (0.2) 
and a certain group D, which commutes with C and D: | 
Dp oe cÓy 31:(0.8) 


where C denotes the charge conjugation operator. 

We define as isotransformations with respect to a given equation these trans- 
formations, which are permitted by this equation and which commute with the proper 
Lorenz group C. From this definition it is seen that (0.3) is the isotransformation 
group connected with equation (0.1). l 

If we introduce in (0.1) the electromagnetic field 4,, we obtain 


POR ie mmy 0 (0.4) 
(499) 
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From the assumption that charge is a scalar-isoscalar quantity it follows that the 


group of isotransformations reduces for (0.4) to 


v" = av (0.5) 


For a = 1 (0.5) is the well-known gauge transformation group. 

In this paper we enlarge the isogroups, which are permitted by equations (0.1) 
and (0.4) 

It is well-known that the 4-spinors admit a second (besides C) 6-parameter 
group €' (Gursey 1958, Rzewuski 1958). This group commutes with C and can be 
simply written as follows (Rzewuski 1958): 


* 
e Com dcr bc; 
Y. i; 
Cee + de; 


ad be] (0.6) 


where c,, c; 


For Dirac's 4-spinors €’ can be expressed in a form independent from the choice 
of the y,„-matrices representation (Tuokoka 1959): 


are Weyl's 2-spinors. 


C: y —(A + By y + (C+ Dyg Cy : (0.7) 
|4|? — |B|? — |C|? + |D? 2 1 AB* = CD* (0.8) 
The group (0.3) restricted by the condition 


la — |t = 1 

is a subgroup of €”. 

The product € x ©’ was investigated by Gursey (1958) and Rzewuski (1958). 
In Gursey's paper this group is applied to the nucleon equation. He shows that the 
masses of the proton and neutron remain scalar-isoscalar. We apply € x C’ to Dirac’s 
equation (0.1) and (0.4). It is easy to show (see Sect. 1 and Sect. 2) that mass m and 
charge e must be replaced by the operators M and E, which transform with respect 
to C’ as an isovector and an antisymmetrical isotensor of second rank. Equations 
(0.4) appear then as the general equation containing M and É written in a special 
isorepresentation of the mass isovector and charge isotensor. 

In the general equation, in the case A, = 0, the mass isovector components x i 
are included symmetrically with the momentum operators id „(see Sect. 6). This sym- 
metry is most clear if we go over to the corresponding second rank equation: 


(O — xi) y = 0 (0.9) 


‚It is seen that the momentum in the usual space and the mass vector components 
in the isospace appear in (0.9) symmetrically. After some restriction, which can be 
simply interpreted, we obtain from (0.9) the generalized Pais equation, introducec 
by Rayski (1956). | 
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The charge conjugation cannot be realized by means of €" (cf. (0.8)). It is obtained 
as a product of two isoinversions. In a special isorepresentation of our equation, cor- 
responding to Dirac’s equation, the first isoaxis and isotime are inversed. 

If we consider the 4-spinors as operators, which satisfy the conventional anti- 


commutation relations, these relations are maintained only under the 3-parametric 
subgroup of (0.7): 


y = ay + dy, Cy jal? + |d — 1 (0.10) 


It is the well-known Pauli group, which describes isorotations in a (1,2,3)-sub- 
space of the full Minkowski’ isospace ©’ (Pauli 1957). 


l. The mass as isovector 


It is well-known that we can treat a 4-spinor as a 2x 2 matrix y (Gursey 1956). 
For Weyl's representation we obtain 
Yaı C, Tau = Cs (1.1) 


Under a proper Lorenz group € this matrix gets multiplied on the left with 
a 2x2 unimodular matrix. We can multiply also on the right by 2x2 matrices. If 
we take these matrices to be unimodular, we obtain a 6-parametric representation C’ 
of Minkowski's isogroup. We indicate the elements of y with an index before thc 
semicolon — which belongs to C, and an index behind the semicolon — which belongs 


to C'. 
The Dicac equation (0.1) can be written with help of (1.1) as follows: 


Veh; (os = Fe M,, (1.2) 


where 


M = fe ij (1.3) 


In Dirac’s equation the mass is a scalar-isoscalar. Thus we can only take into 
account these transformations, which do not change the matrix ,M. From (1.3) we 
see, that Dirac equation (1.2) allows only the following subgroup of C”: 


: a b 
Cia) E i a (1.4) 
If we take into account the full group C’, we must treat M as some isospinor 
of second rank with upper indices. Thus we can generalize (1.2) as follows: 
ab; NER = aT, gra (1.5) 


where a — arbitrary complex scalar-isoscalar, 
M — arbitrary 2x2 hermitean matrix. 
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After choosing 
ab; — (0 29,5; p ase In; (1.6) 


Unc un 
M. 


iB T (0; i d ? 50) sap = Oi Ze (1.7) 


where o,., o, are Pauli matrices, x, o. — real, we obtain that 
> 3 


Ms MR} — (chy — x) 8,07 | (1.8) 
and 
«59, = — 10%, (1.9) 


From (1.5) and its complex conjugate 
gef, = a* (esty M. (1.10) 
we obtain the second order equation 
O Fas = lal? (d — 2%) eG (1.11) 
Comparison with Klein-Gordon equation yields: 
la]? (22; — 22) E (1.12) 
In further considerations we put 
le? (1.13) 
and we have 
Wa Hq cm (1.14) 


It is seen that 'P,., from (1.5) or (1.10) satisfies also the Klein-Gordon equation 
with the mass m if Ms4v describes a space-like mass isovector %.„» the length of which 
is determined by (1.14) and equals m. 

The usual Dirac equation (0.1) or (1.2) follows for 


xi = (0, 0, — m, 0) (1.15) 


as may be seen from (1.3) and (1.7). 


It is easy to show that (1.4) or (0.3) generates rotations in the (1,2,4)-pseudoeucli- 
dean subspace of full isospace. 


2. The charge as isotensor 


It is well-known that in presence of the electromagnetic field we must introduce 
the following change in Dirac’s equation (see (0.1) > (0.4)): 


Qu; a Ou; — ie Ap | (2.1) 
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Because C’ transformations (0.7) do not commute with ie, we see that following 
change 


Q ap; > ap; + € Aap, (2.2) 


where c is scalar-isoscalar and 


Ash; = Gab; Aus (2.3) 
is not equivalent to (2.1). 
A second possibility is 
Pai; > Paps + Aap; Fd (2.4) 
We can write (2.4) also in the form: 
Pab; = Pab; 3 A ab; Ld (2.5) 


Equation (1.5) has now the form invariant with respect to the full C x C” group 
and includes the electromagnetic field 4,;.: 


Des: (uv) ds Agi, (659 p =F) DE Me (2.6) 
We can express Eg” with help of an antisymmetrical tensor e,,”: 
Ej. TE eu O50. 0,57 (2.7) 
If we express (2.1), which acts on the 4-spinor y, with help of the 2X 2 matrices, 
we obtain: 


ENS E (2.8) 
where 
e 
0 — og 00 
(e;2,)) = E 000 (2.9) 
0 000 
0 000 
Assuming 
Cute =e (2.10) 


we obtain (2.8) as a special representation of the general (2.7). It is Dirac's represen- 


tation of the isotensor e,,,,. 
From (2.7) we see that 


Lo = (ie,; um ai d.i o)? à (2.11) 
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i in (6,25; — €4us € 13) (2.12) 
i d; = (6510: €:99* € 39) 


Condition (2.10) can be written as 
e? — d$ =e (2.13) 


Bo 
For the unitary subgroup (0.10) of C’ c'* has the same trans formation properties 
as c,,. Thus for rotations in (1,2,3)-isosubspace e,, is a pseudoisovector and d,;— 


an isovector. For Dirac’s case we have: 
eD—(0,0,—e) d? =0 (2.14) 


We see that in the isospace corresponding to C’ we have a similar situation as 
in usual space. If we write the following analogy: 


where E; and H., are components of electric and magnetic strenght vectors, we 
obtain the following correspondence: 


p cof rise Susi, A | (2.16) 


uv 
These analogies can be useful dod if we want to write in € s some equations of 
motion. à 


3 Charge conjugation 


It is inc to mention that charge conjugation is not ied in (0.7). 


Charge conjugation C can be expressed in Weyl's 2-spinors: 


£s EC; Zo - xs 
C cem am EM .: 7. 


a 3 H trali 


From (1.1) we see that 


(OM OR 
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4. Description of equations (1.5) and (2.6) in 4-spinor notation 


Cs EI 
Yr (3) = (ve (4.1) 


Put 


2 Pa; 
(y= y= pas (4.2) 
Using (4.1), (4.2), (1.6) and (1.7), we can write (1.5) in the following form: 
(y, 9, + M)y =0 (4.3) 
where 
M = Msg — 9:9 Ys — (651 + 19,9) ys G (4.4) 
and 
0 —a*o, 0 —ia* 
ULM e 0; 3 ovs E 0 ) (45) 


Thus instead of the usual Dirac equation we obtain an analogous equation, 
where M is an isovector (4.4) invariant with respect to Lorentz transformations (0.2). 

Equation (4.3) can be written in an another form if we introduce the following 
8-component spinor ®: 


ó— | Je) (4.6) 
In terms of ® equation (4.3) takes the form 
(,8,- 1x) Ð =x 50 (4.7) 
where 
D,—y,X I9 s —I? x a, (4.8) 


I? — unit matrix of n-th rank, 
X  — 3-component mass-isovector in (1,2,3)-subspace of C'. 
Now we introduce the general charge isotensor (2.7). Calculating with help of 
(4.1) and (4.2) expression (2.7) we get instead of (2.1) 
9,>9, + EA, (4.9) 


where 


2 È= le 12 + ys Cen +19 Gea an Cei Siy i Cg. + Ys €,34 (4.10) 
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If we insert in (4.10) Dirac’s isorepresentation (2.9) of e,,,, we obtain the usual 


formula (2.1). 
By virtue of (1.6), (1.7), (4. U, (4.2) and (4.5) equation, (2.6) can be written in 


the following form: 


(y, (0, + EA) + M) y —0 (4.11) 
where É and M are determined by (4.4) and (4.10) with additional conditions a. 14) 


and (2.10). l 
Equation (4.11) may be written with help of the Serine ®. After simple calcu- 
lations we obtain 


T,1,+i@—il,dt-Al tr BO=xy 6 (412) 
with 
ee E 
M; =V% — m? (4.13) 


In Dirac’s case we have: 


(D. (0, > ders A) ma) 9? 0 = (4.14) 


rc T. 


De GREEN of (4.12) to Dirac's equation. 


Let us write (4. 12) in a particularly simple isorepresentation, From (4. 13) it is 


seen that (6, d) is a six-isovector with positive invariant form. Because ® transform — 2 
with respect to C’ according to ; 


Bean oi freal See ae 
& and f can be chosen in a such way that —— eee 2a 


QAI. m CG ‘Gd o agi 
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This pseudoeuclidean rotation transform (5.3) into 


(59, — en 4) T KY} © = 0 (5.7) 
where 
R = Ku 
X. = x (5.8) 


After a second rotation 


.« 7 
ut AX. 
D” =e 2° p" iga= (5.9) 
2 


which does not change the charge operator, we obtain finally: 
Xa =0 x4 —Ywg dox Xx. Hs (5.10) 
Equation (4.12) can now be written in the following most simple form: 
U, (8, — te t3 A,) + Ta 4.5 + 19 4) Ó = 0 (5.11) 


It is seen, that equation (5.11) is essentially different from (4.14); In Dirac’s 
case x, == 0. It is seen further that in the general case the mass-isovector % is not 
parallel to the charge-pseudoisovector e. 


6. Symmetry properties of equation (4.7) 
If we take in (4.8) the Weyl representation of the y -matrices,, 
Vi= 602; Va Or (6.1) 


we obtain from (0.2) the following transformation law of the 8-spinor ® under C-roia- 
tions: 


D' — eA cire & E (6.2) 
Here A, B — two arbitrary real 3-vectors, and 
0 —35 x I9 (6.3) 
Comparing (6.2) with (5.1) we see that transformation K: 
P —> P = } (1 + 0?) D = KÖ (6.4) 


has the following property 
: KCOK1-—C' 
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and 
KO, Ke C (6.5) 
because of K ! = K. 
With help of (6.1) equation (4.7) can be writen as follows: 
fo, (2A — yd) x 19 + tx Iruo} 6 = 0 (6.6) 
Performing here transformation (6.3) we obtain 
(o, x 10 - Ep — I's) + (Ex Yon) x I9) 6 —0 (6.7) 


or, introducing the /', — matrices 


(P, s, + Tp) P = 1,9, 0 (6.8) 


It is seen that in (4.7) the roles of the matrices J", and (z;, ‘I’;) are interchanged. 
Formally we obtain (6.8) from (4.7) through the following substitution: 


EEE (6.9) 


This is connected with the fact that the choice of C as space-time rotation group 
and C’ as isospace rotation group can be reversed. 

If Ø describes plane waves with momentum p,, in usual space-time, the cor- 
respondence (6.9) may be written 


Yay = ipy; (6.10) 


Equation (4.7) or (6.8) describes plane isowaves in (1, 2, 3) — isosubspace. The 
isodirection and isovelocity of these waves are determined by the isomomentum x, 
Because x,, is isospatial (see (1.14)), the isovelocity must be smaller as on the iso- 
-lightcone. For example, Dirac’s equation is an equation for standing isowaves directed 
paralell to the third isoaxis. 


7. Charge gauge group 


Let us investigate the general equation (4.12) written in the isorepresentation 
(5.11). It is seen that the following two l-parameter groups of isorotations do not 
change the charge operator 


G: Go D = en p 
G,: D> D = Pp (7.1) 


G, describes rotations in the (1,2)-isoplane, G, — pseudoeuclidean rotations in 
the (3,4)-isoplane. 

The group of transformations, which do not change the mass parameters, is 
given by the 3-parameter group M of rotations in the 3-dimensional isosubspace, 
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perpendicular to the mass isovector x (0, x.o, %.3) from (5.11). We have 
M: Ó-—d' eo To n) iT, 9 (7.2) 


where T, indicates the new basis: 


l l 
T; — E Be (X55 T — 2a t3), x: (Kar, + 2 (7.3) 


T, and T; are orthogonal to the mass vector, T, is to this vector parallel. Generally 
we can express this new basis with help of the following expressions: 


ITA RR 
Ty = -axa xj- 
T, = a)r (7.4) 


LE YR T 


where a is an arbitrary unit 3-vector, not parallel with x. To obtain the mass- 
-isovector occuring in (5.11) we put a =) l, 0, 0). 

The charge gauge transformation group is a continuous isogroup which does 
not change the mass and charge operators. It is given, therefore, by the common 
part of G, and M. It is easily seen, that only in the two cases x, = 0 or x,y =0 we 
can define continuous charge gauge transformations. 

In the case x, =0 (Dirac's case) we have 


WARES (7.5) 
and from (7.1) we see that G, =M if 
Deu 0 “=e (7.6) 


G, is the conventional gauge transformation (0.5). 
In the second case, if 


(0,7. 070) (7.7) 
we have 
T, = 45 
Taking 
qe oc, = Oke (7.8) 


we obtain G, = M. 

In the general case we cannot define the 1-parameter gauge transformations. 
Thus, from the point of view that a continuous gauge group is obligatory, we obtain 
only a second possibility beyond the Dirac case. This second equation can be written 


510 J. Lukierski 


in terms of the 4-spinor y as follows: 
Yu (9, — ieA,) v = imysy' (7.9) 
The gauge transformation group is now 
y= el y (7.10) 


It is seen that equation (7.9) is under (7.10) invariant. 

Finally we see, that mass-isovector x and charge isopseudovector e must be 
parallel (eq. (0.3)) or perpendicular (eq. (7.9)). 

Because in (5.11) the mass operator is invariant under isoinversion laws 
also the charge conjugation operator for the general equation (5.11). Simultaneous 
mass and charge conjugation correspond to the following isorotation: 


®' = it,® (7.11) 


This is an rotation about and angle zx in the (2,3)-isoplane (the mass-vector is 
situated on this plane). Thus mass conjugation is realised by help of the total reflexion 
Th34 of the isospace: 


lene (Uy CPU (7.12) 


By virtue of (4.12) we obtain the mass conjugation as strong reflexion in isospace 
or in usual space-time. 

We can restrict the 6-parameter isogroup C’ to Pauli's group (0.10). Because 
only this 3-parametric subgroup of C’ preserves the conventional anticommutation 
relations (Pauli 1957) we must, in quantum field theory, either change these relations _ 
or restrict the group (0.7) to (0.10). 

If we put in (4.7) 


4.9 =0 (7.13) 


we obtain commutability of isotransformations with the full Lorentz group (including 
inversions). This is connected with the fact that inversions in space-time cannot be 
realised without iso-inversions, for example isotime inversion. If we demand isotrans- 
formations to commute with the full Lorentz group, we obtain a second reason to | 
restrict C' to Pauli's group. A restriction of this type excludes G, and admits only 
Gi, gauge transformations. Dirac’s equation written in the Heg isorepresentation is, 
therefore, the only possibility in this case. 


8. The connection with non-local theories 
Equation (4.7) can be written in a symmetrical way, if we introduce a new coordi- - 


nate-space 7,,, corresponding to the isomomentum space 


5 | 
X su > i EX 
os (8.1) 
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From (4.7) one obtains in this way the following Schródinger equation 


.ID 
I el (8.2) 
where 
dod AE’ Sr 
(sr m 2 A X; 
AEN OX 4 A=4, 5, 6,7: X e 7, 
and 
ale 207r (8.3) 
We have explicitly from (4.7): 
9 . d 9 
Zi d Ei (8.4) 
The second order equation, corresponding to (8.2), is 
(ERS, 0 (0) — 0 (8.5) 


The two spaces x, and n, are independent. 
If we take into account the arguments of Section 7, we can restrict the four- 
vectors x,, N, to their parts in (1, 2, 3)-isosubspace. From (8.5) we obtain that 


(pu + x2) Ø (x, N; p, x) = 0 l (8.6) 


This is the equation proposed in Rayski’s generalization of Pais theory (see for- 
mula (21), Rayski 1956). _ 

The author is indebted to prof. Rzewuski, prof. Rayski and prof. Lopuszanski 
for helpfull and valuable discussions. 
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QUENCHING OF PHOTOLUMINESCENCE OF SOLUTIONS 
BY NON-ABSORBING FOREING MOLECULES 


By Jan GŁOWACKI 
Physics Department, Normal School, Gdansk 
(Received January 22, 1960) 


The theory of luminescence quenching proposed by A. Jablonski in 1954, has been 
checked experimentally for liquid and rigid solutions of fluorescein, rhodamine B and 
*yellowish* eosin quenched by KJ and tartaric acid. This theory with a simplified model 
of the centre of luminescence agrees very well with the experimental results of the meas- 
urements of photoluminescence quenching. It has been also shown that, if the molecules 
are ions, the probability of transition of a quencher from one group to another depends, 
besides the viscosity of the solution and its temperature, on the charges of the luminescent 
and quenching molecules and on the dielectric constant of the solution. 


8 I. Introduction 


In the theory of quenching of photoluminescence of solutions by Jablonski (1954) 
the centre is considered to consist of an excited.luminescent molecule enveloped by 
shells consisting of monomolecular layers of the solvent, which may contain quench- 
ing molecules. The quenchers may change the shells. The probability of a transition 
of a quencher from one shell to another depends on the viscosity of the solution and 
on its temperature (Jabłoński 1956). If the solution is rigid, the quenchers remain _ 
infinitely long in their shells. The probability of quenching of an excited molecule 
by a quencher depends on the shell in which the quencher is “situated. 

The mean value (i. e., the observed value) of the relative yield of luminescence 
for solutions containing quenchers was calculated on the basis of a simplified model 
of a luminescence centre (Jabtoriski 1954, 1957, 1958). Assuming that the quenching 
probability is different from zero only for quenching by quenchers situated within 
a certain active sphere, the yield for a centre containing k quenchers (mutual diffusion 
of quenchers and luminescent molecules neglected) is (provided the volume of the 
sphere is suitably chosen) 


ee a) 
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where n, denotes the yield for k = 0. By applying for k the Smoluchowski distribu- 
tion (1917) 


k 
Pise. : (2) 


DN 
ki? 


the following expression was obtained 


foe) re» 
TV Bu-- (3) 
No er! 


v 


where » = nv (n — nümber of quenching molecules per unit volume, v — volume 
of the active sphere). 

The aim of the present paper is to check experimentally the above simplified 
version of the theory for a case of quenching by non-absorbing foreign molecules. 


§ 2. Quenching of photoluminescence of rigid solutions 


In order to check formula (3) solutions of dyes were used. Rigid solutions of 
fluorescent substances were obtained by using sugar syrup of great viscosity in which 
fluorescein and KJ (quenching substance) were dissolved. To estimate the order of 
magnitude of the diffusion coefficient in such a viscous solution, the simple, approxi- 
mate Herzog-Polotzky formula (1914) 


D = ya 86400’ (4) 


was used, where u denotes the molecular weight. The very accurate experimental 
results obtained by Hodges and La Mer (1948) for fluorescein in water (D = 3.7 x 
x 10-6 cm?/sec) agree well with those calculated by means of the Herzog-Polotzky 
expression (D = 4 x 10-6 cm?/sec; for other solvents D is inversely proportional 
to the viscosity). Assuming the mean life of the excited state of the luminescent 
molecule of fluorescein to be v ~ 5 x 10^? sec, we obtain for the mean diffusion 
displacement during v the value 0.1 A; the radius of the solvent molecule is r ~ 5A. 
Therefore, the solution of fluorescein may be considered with good approximation 
as a rigid solution. The results of the measurements are given in Fig: 1. It may be 
seen that expression (3) fits very well to the experimental results. 

The most interesting seem to be the results of the measurements of the fluores- 
cence yield of rhodamine B in ,,sugar glass" quenched by KJ. It is worth noting that 
according to Vavilov and Levshin (1926) and Pringsheim (1949, p. 448) the rhodami- 
nes seem to be the only exception, being non-phosphorescent in rigid solutions at 
room temperature and emitting no slow fluorescence at the temperature of liquid air. - 

Rigid solutions of this dye were prepared by melting sugar with a small quantity 
of water and subsequently pouring this syrup into a metallic ring placed upon a glass 
plate. After several minutes, when the solution becomes very viscous, the ring were 
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covered with a glass plate. In this manner almost rigid solutions were obtained. These 
solutions were investigated at a thickness of the layer amounting to 0.3 cm. The 
results of the measurements are given in Fig. 1. 

In all cases the photoe ctric method with photomultiplier was used for measur- 
ing the relative yield of fluorescence. The solutions of fluorescein were excited through 
a blue filter and the solutions of rhodamine B — through a green filter. Fluorescence 


Te 


g% 7 2 v 


Fig. 1. Relative fluorescence yield versus v. 


— theoretical curve, 
©  — expermental points for fluorescein in sugar solution quenched by KJ, 
€  — experimental points for rhodanine B in sugar solution quenched by KJ. 


01 02 0.3 0.4 05 
Concentration in g/cc 


Fig. 2. Relative fluorescence yield as the function of concenctration of quenchers (KJ). 


©  — experimental points for fluorescein in sugar solution quenched by KJ, 
@  -— experimental points for rhodamine B in sugar solution quenched by KJ, 
—— — — theoretical curve ‘calculated from Jablonski’s theory, 
— — — — theoretical curve calculated from Frank-Vavilov's theory. 
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intensity for fluorescein was measured through a filter in the spectral range of 5100. 
to 5350A and for rhodamine B — through a filter in the spectral range of 5600 to 
6300 A. 

It will be seen that the simplified model of the luminescent centre suffices for 
the description of the quenching process of fluorescence by non-absorbing molecules. 
In Fig. 1 the experimental values and theoretical results corresponding to R = 5.74 A 
for fluorescein and R = 6.8A for rhodamine B are compared. The results of the 
measurements show that quenching is caused, most probably, by quenchers situated 
in the shells nearest to the luminescent molecule. Unfortunately, it was not possible 
to exceed the concentration of 0.5 g/cc of the quenchers. 

The experimental results of quenching of rigid solutions are compared also with 


Frank Vavilov’s theory (1931)!. The results of measurements show that, most probably, 


the quenching does not proceed according to a simple exponential law (ce, Fig. 2). 


§ 3. Quenching of photoluminescence of liquid solutions 


Although the mutual diffusion of quenchers and luminescent molecules and the 
orientations of their molecular axes are neglected in expression (3), this expression 
agrees also with the experimental results of quenching by absorbing molecules of lumi- 
nescence of non-viscous solutions. In order to check expression (3) for the quenching 
of luminescence of liquid solutions by non-absorbing foreign molecules, aqueous 
solutions of fluorescein quenched by KJ were used. The results of the measurements 
are given in Fig. 3. Expression (3) with the radius of the active sphere R = 17 A 
fits to the experimental results very well. 


Wn, | 
7-9 


0 7 2 3 4 5 6 


E 
Fig. 3. Relative fluorescence yield for fluorescein in water quenched by KJ versus v. 


— theoretical curve, 
©  — experimental points. 


* + : Ä 
1) Their formula for the relative fluorescence yield is n/N = exp (—wc)/(L+ kc). If k = 0 (the case 
of rigid solution), 7/7) = exp (—wc) where c is the concentration of quenchers and w a constant. 


Quenching of Photoluminescence of Solutions 517 


The author investigated also the dependence of the degree of polarization of 
fluorescence on the concentration of the quenching molecules. This dependence is 
described adequately by the. Perrin (1926) expression 


Ji il v [3 l 
p Doo To \Po Po 


where p is the degree of polarization for a given concentration of quenching mole- 


cules; p,, — the degree of polarization for total quenching of luminescence; pg — 
the degree of polarization in the absence of quenching; t — the mean life-time of 
the excited state of the luminescent molecule for a given concentration of quenchers; 
To — the mean life-time of the excited state of the luminescent molecule in the ab- 
sence of quenching. 

The well known formula given by Perrin (1929) (valid in some cases only) is 


= — , 6) 
Io No To ( 


From (5) and (6) we obtain the linear dependence of the reciprocal degree of polari- 
zation on the yield of luminescence 


vnl zr (7) 


where 


0 7 2 o 4 9 CN 


Fig. 4. Reciprocal of the degree of polarization of the fluorescence and relative fluorescence yield o 
rhodamine B in water quenched by KJ às the function of v. 
— theoretical curve, 
© and € — experimental points of Frank-Vavilov for relative fluorescence yield and reciprocal 
of the degree of polarization of fluorescence, 
© and ® — author's experimental points for relative fluorescence yield and reciprocal of the degree 
of polarization of fluorescence. 
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The check of the dependence of p,/p as function of v is given in Fig. 4. Expres- 
sion (7) was compared with the experimental results obtained by Frank and Vavilov 
(1931) and by the author at high concentrations of the quenchers for rhodamine B 
in water, quenched by KJ. The author’s results of measurements of the degre of pola 
rization of luminescence light were made by the visual method of a Savart plate and 
Arago compensator. The constant po/p,, was determined graphically from Graph 5. 


0 05 7 9/5, 


Fig. 5. po/p versus m/m. 


8 4. Quenching of photoluminescence 


of water-soluble ,,yellowish” eosin by tartaric acid 


Another aim of the paper was to examine the quenching type of luminescence of 
yellowish” eosin (microscopic dye, molecular weight 687.026) in water by tartaric 
acid. 

Since the tartaric acid is optically active (it rotates the plane of polarization), 
which made measurements of the degree of polarization difficult, such investigation 
has not been made in this paper. The luminescence of eosin was excited with light 
which had pasged through a Zeiss mirror monochromator giving the exciting light 
in the spectral region Ana, ~~4750 A. The fluorescence intensity was measured by 
means of a photo-multiplier (type FEU-19M) through a yellow-green filter of 5300 
to 5900A, Amas = 5600 Å. | 

The results of measurements of the intensity of luminescence and pH for different 
concentrations of tartaric acid are given in Figs. 6 and 7. Fig. 8 shows the absorption 
spectrum. The change of the absorption spectrum with the concentration of tartaric 
acid shows that the quenching of photoluminescence of eosin is the result of the for- 
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Fig. 6. Relative fluorescence intensity of eosin in water quenched by tartaric acid. 
— tehoretical curve calculated from Eq. (8) 
— —— — theoretical curve calculated from Eq. (9) 


©  — experimental points. 
PH 
70 
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Fig. 7. pH of eosin solution in water versus the concentration of tartaric acid. 


mation of non-fluorescent complexes of eosin and tartaric acid molecules (,,static” 
quenching). The experimental points are in good agreement with the theoretical 
curve of hyperbolic type (Fig. 6, curve a) 


" 1 
Dd Ke) 9) 


The insignificant systematic divergences between the experimental points and the 
theoretical curve at high concentration of the quenchers (tartaric acid) are probably 
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due to the effect of the fluorescence of complexes. Taking the fluorescence of com- 
plexes into account (Férster 1951, p. 203), we obtain a new expression 


"A SI Kee 
4500 5000 5500 6000. A 
Fig. 8. Absorption spectrum of eosin in water. 
1 — solution wtih c = 8.4x 10-8 g/cc concentration of eosin without quenching molecules (tartaric 
acid). 
2 — solution of eosin with c = 0.472 10-? g/cc concentration of tartaric acid. 
3 — solution of eosin with c = 1.148x 10-? g/cc concentration of tartaric acid. 


4 — solution of eosin with c = 7.84 10-? g/cc concentration of tartaric acid. 


where A and K are certain constants. Expression (9) agrees better with the experi- 
mental results (Fig. 6, curve b). 


Another interesting feature of the quenching of eosin by tartaric acid, studied in 


this paper, is a linear dependence of intensity of luminescence on the pH (Fig. 9). 
This dependence is given by the following empiric expression | 


H s 
—- +C= BA (10) 
where B and C are constants, Jand Jọ denote the intensity of luminescence in the pres- 


ence and absence of quenching, pH and pH, are the pH of the solution in the pres- 
ence and absence of quenching respectively. ` 
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0 0.5 | 7 pH/ph, 


Fig. 9. I/I, versus pH/pH,. 


8 5. A note on the constant A,, for the shell model of the luminescent centre 


If the transition from one group of centres to another is taken into account, 
the behaviour of the system is described by a system of interdependent differential 
equations (Jabłoński 1954), in which case the constant 4; denotes the „transition 
probability” of a transition of a centre from group i to group j°. If we assume the 
simplest case in which all the centres are divided in two groups only, the centres 
with one single quencher in their shells forming group 1, and those with no quench- 
ers — group 2, then A, is the probability of a transition of a centre from group 1 
to group 2, and A, — the probability of a transition of a centre from group 2 to group 1. 
A, dt denotes, therefore, the probability that a quencher enters into the active sphere 
of a given centre from outside during dt. It is, therefore, the process of diffusion 
which was investigated by Smoluchowski (1917). According to Jabłoński (1956) the 
following expresion was obtained 


dy = 4a DRN i (11) 


` by assuming that 2 RV «JV aD <t (t is the mean life of the excited state of the lumi- 


nescent molecule), where 


er 


6an \r r 


2) All the centres are divided into groups according to the occupancy of their shells by quenching ` 
molecules: the contres with one single quenching molecule in the first shell belong to groupl, the centres - 
“with two quenching molecules belong to group 2, and so on. i 
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Dee eee 
is the diffusion coefficient of the luminescent molecule and of the quencher; R is 
the radius of the active sphere; N — the concentration of the quenchers. Hence 


(Jabłoński 1954) 
jM bare (12) 


where v is the volume of the active sphere and v = Nv. The expression A, depends, 
therefore, on the viscosity 7 of the solution, on its temperature T and on the radii 
of the luminescent molecule and the quencher. The expression (11) can be correct 
only if the fluorescent molecules carry no electric charges. If the molecules are ions, 
as is very frequently the case, A, depends on the Coulomb electrostatic force. The 
diffusion problem in this case was investigated by Smoluchowski (1915), Chandrasek- 
har (1943) and Debye (1942). According to Debye, if the molecules are ions, 4a DRN 
must be multiplied by 


6 
= l 
TENE LEN (13) 
where 
1 Ze 
ô = TR (14) 


Here, Z,e and Z x denote the charges of luminescent molecules and of quenchers, 
e — the dielectric constant of the solution. If Z, Z, e? — 0 (one or both molecules are 
neutral), then A = 1. 

By taking into account eq. (13), we obtain instead of (11) 


À, = AxDRNA. (15) 


Expression (15) depends, therefore, on the charges of molecules and on the dielectric 
constant. 


8 6. Concluding remarks 


From the examples investigated in this paper the following conclusions can be 
drawn. For an adequate description of luminescence quenching the simplified model 
of luminescence centre with the active sphere is sufficient. The experimental results 
of photoluminescence quenching by non-absorbing foreign molecules in liquid and 
rigid slsutions agree with Jablonski’s theory. It describes also adequately quenching 
by absorbing foreign molecules the absorption spectrum of which overlaps the fluores- 
cence band of the luminescence molecule (Jabłoński 1957, 1958, Głowacki and Po- 
hoski 1959). 

The author is greatly indebted to Professor A. Jabłoński for suggesting the 


subject and his numerous friendly criticisms. He also wishes to thank Dr. J. Szych- - 


linski, Mr. T. Latowski, and Mr. Z. Polacki for their assistance with the measurements. 
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The angular asymmetry of neutrons from the capture of polarized w-mesons in 
complex nuclei is estimated, assuming the applicability of the Fermi gas model of the 
nucleus. It is shown that the four-fermion V—A interaction.together with the conserved 
current hypothesis lead to considerable asymmetry. 


l. The idea of universality of weak interactions leading to processes with four 

fermions expressed ten years ago[!] has got in the last two years strong experimental 
and theoretical support. The notion of the universality was made more precise. Namely, 
it does not mean that these processes may he described by a set of couplings with 
the same constants C s» Cy, ete.; the universality is expected for a primary weak 
interaction, which may lead to different effective four-fermion couplings in different 
elementary processes as a result of an internal structure of strongly interacting particles 
entering in these processes. It is known that the experiments on f-and w-decays 
are consistent with V-A coupling. To explain non-renormalization of the vector 
.coupling constant in ß-decay, Feynman and Gell-Mann[?] supposed that into the 
vector coupling with (ev)-curent enter not simply the nucleon current, but a total 
conserved isovector current. This hypothesis leads to additional physical consequences: 
it gives corrections to the energy spectrum of f-electrons[?], introduces the angular 
correlation B—y in allowed transitions[4], etc. These effects are, however, relatively 
small because of the small nucleon momentum transfer involved in f-decay (g? = 
= (p, — py? = 0)". 

If the idea of the universality is correct, then the least studied up to now x~ -capture 
proces must be described by the V-A interaction with the conserved vector current. The 
effect of the latter — if it exists — must manifest itself here more significantly because of 
the larger momentum transfer (q? ~ 0.9 m?, m—the muon mass). In the present 
paper the process of muon capture is considered from the point of view of the exami- 


nation of V-A coupling and especially of the conserved vector current hypothesis. 
a a ee a eR 
*) We adopt h =c =1 : 

(525) 


526 W. Majewski 


2. We assume 


Gye G A ; 
Hise = A Ji [Prya (l + v9 Yal — Vi jt [Pri ya ys d. +75) vu] + hic., 


A = vae yu t ib T Va Pn — (s q,) T- Pa] + (other posisble terms), 
Ji = Yn TU, Ys YN- 


The vector current J“ fulfils the condition iy ; 0. To the first order in the weak 
interaction the S-matrix element for muon capture by a proton has the following 
form 


See er) 04 (peepee pyar, R: 
oye <n jilp [u, ya (1 ys) uel — = <nlalp> [u,i ya ys (l+ yg) ul- 
y2 y2 


Here |n>,|p> denote the states of physical nucleons. The effect of strong interaction 
is taken into account phenomenologically by writing the nucleon matrix elements in 
a general form following from the Lorentz-covariance. Taking into account the charge 
symmetry and the T-invariance of the strong interactions, we have only four form- 
factors[5]: 


x 


<n|jilp> = fy Un |» "t 2M Ore (Pn - poe Up 


A = : a 
<aljalp> = fy Un E i dtd gulam (ps — pois | Up; 
M is the nucleon mass; u,,,, are bispinors of the corresponding particles, 
‚which are for the present assumed to be free; fp, x, A, & are functions of q?. Using 


the conservation law of the 4-momentum and the Dirac equations for leptons, we 
obtain l 


G 
m = fel ta (ri app me P») e | n vd + 


+A [u, Va ys uj] [u, Ya (1 + 75) Yul + € [u, Ys us] [u, (1 — ys) 2 : 


We see that taking into account of the effect of strong interactions changes the effective 


V and A interaction strengths: G/V2 > fy G|V2, — G|V2  — å fy G/y2, and intro- 


duces two new couplings: a pseudoscalar one with a factor cS? — —£ fy CIV2, and 
a derivative coupling — the so called weak magnetism term — proportional to 
f x G 
ENN ty —. 
2M *" y2 


w- wn 
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The conservation of the vector current j? and the charge invariance of the strong 
interaction lead to the equality 


<nljilp> = <N|ji|N> = 
= um |». (Ep Fo) + 514 (n^ Fp — wt Fg) Sip (pa — pe u (py). 


Jj is the isovector part of the charge-current density of nucleons and pions; 
x? = 1.79, x" = — 1.91 are the anomalous magnetic moments of proton and neutron, 
respectively, in units e/2M; 
Fp", FP" are the nucleon electromagentic formfactors, known from experi- 
ment [6] to be 
ERES LE pes Felt <r 2g? — 0,7 


for the momentum transfer involved in muon capture, if the root mean square radius 
of the nucleon. charge density is <r> = 0.80x10-!3 em. The formfactors in 


p-decay and w-capture were investigated in[?] by the method of dispersion 
relations; the result is 


fy Mop: ~ fy A ao and (fy &)/( fy A) | 2=0.9m: ~ 8. 


The formfactors fy, etc. for q? = 0 are measured in f-decay, therefore, using the 
experimental values f (0) &1, 4(0)~1.25**), one obtains for muon capture. f,= 0.97, 
x = xP—X" = 3.7, A = 1.29, Æ = 10.3. Let us notice that if we do not assume the 
conservation of the vector current, then 32«0.25[?], and the weak magnetism term 
is negligibly small in muon capture too. 

Rejecting the terms ~ (v, 4/c)? ard assuming the muon to be at rest at the moment 
of capture, we obtain 


G l Wet ES > | ;5" > 
Bis ya" (ent) (a+ b «8 x) 


acl P b. cL | 
ges DW ees + Pot Pn 4g Po— Pn A N N (1) 
m rh SV 2M 2M 


y=p MB EZ o> aes o> p, are matrix elements taken between the spins states 
of nucleons and leptons, respectively; u = 1 + x? — x". 

3. In [7] it was shown that the effect of the conserved vector current (which redu- 
ces essentially to an about fifteenfold increase of the coefficient x) considerably changes 
thé angular asymmetry of the neutrons from capture of polarized muons in hydrogen. 
With inclusion of an effective pseudoscalar term (omitted in [7]) we obtain from 
(1) for capture in hydrogen me = vm 


e Baits es alt aÈ, ñ) ph do, 
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I —1--3224- B(L + 22 + 24p.— AE) -- 1 B2 (1+ a2 + 292 + E2244) (2) 
Ig = 1— 42 + B(1 + a2 — 2d — AE) -- 1B3 (1 + APH 2? + £2— 228), 


> 


SE (here and further on p, =. p= = j P, is the polarization vector of 
Pn 


captured muons). The space part of the muon. wavefunction is taken as a solution 
for the K-orbit at zero: |o(0)|? = l/mag, ag = (me?) . Numerical estimation 
leads to the result that the effect of a conserved vector current increases the 
capture rate by 17%; still greater is its influence on the asymmetry parameter a: 
it increases by 33%, from —0.33 to —0.44. (This & is calculated neglecting the 
hyperfine structure of the u-mesic atom). 

As pointed out in [8], the experiment on the neutron asymmetry is practically 
unfeasible in hydrogen — there occurs a total depolarization of muons with transi- 
tion of u-mesic atoms into the singlet state from which the capture goes on. It follows 
from (1) that 


. Using] ~ SUN: 
J2 2404-1) ABI 42 (9-1) — 222 E4248] 4 
tiP iW + 2h +]. 


It follows that the conserved vector current hypothesis. increases the singlet capture 
rate by 11% in comparison with the usual V-A theory. 

4. At present only experiments on u”-meson capture in nuclei heavier than Be are 
being performed. For this reason it seems worthwhile to estimate the expec ted 
neutron asymmetry for capture incomplex nuclei: we are interested in the investigation 
of a decrease of this asymmetry in comparison with the hydrogen case. It may be 
emphasized that the principal contribution to & becomes from the relativistic terms, 
namely the pseudoscalar and weak magnetism ones (in the nonrelativistic limit one 
obtains from (2): @ = —0.11). We adopt for the nuclear proton the effective inter- 
action hamiltonian yielding for the free proton the matrix element (1). Our estimation 
of the asymmetry will be probably only qualitative: the proton motion inside the 
nucleus is regarded according to Fermi gas model which already was applied previously 
in-muon capture calculations (see e.g. [9, 10], where the validity of this model is also 
discussed in detail). . 

Namely we suppose that nucleons move independently as free particles in the 
sphere of radius R, being in potential well of the depth Vj = Ep + Ep, and occupy 
all the momentum states up to Pr determined from the condition 


**) See e. g. J. A. Smorodynski Uspekhi Fiz. Nauk, 68, 653 (1959). 
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(for the sake of simplicity one takes Z = N = A/2, Z and N are numbers of protons 
and neutrons in the nucleus, respectively, with A = Z + N); R = R, A^, Ry = 
= 1.2x 103 em, E, = P?/2M ~ 33 MeV, Eg = 8 MeV — the binding energy 
of last neutron. The calculation consists in an integration over the momentum 
distribution of protons of the probability 


dp, [4 LA 
dw = 22|9R[ ETE xP) 5 Ppo/dE. 


It may be shown that the energy — momentum conservation and the Pauli principle 
for neutrons impose on the proton momentum — for given neutron jkinetic energy 
inside the nucleus E, — the lower limit 


Py > Pmin = (M + V2M. Wc oa m 


Because of the approximate character of the applied model only principal relativistic 
terms were retained in the calculation ( that is, terms with the coefficients u and £). 
For the nucleus of spin 0 the following angular distribution of neutrons with respect 
to the muon polarization was obtained 


dw ~I (E,) [1 + « (E,) P, 3] dE, dQ, 


I (E), = jae229 5 ren £) all 


M 1B R 
= — — 2 —— —— — — 
I (Ej) a (Bn) = — |o 22) | ay 0x4 m) 


l deg TS £2 
cope Sedes sp 


kı 


^ kh = M + V2ME,— [(M + V2ME,? — 2 Mm]* 


ky = Eg--m— E 


It follows from kinematics that the emitted neutrons have kinetic energies outside 
the nucleus up to about 16 MeV. I(E,) is the neutron spectrum, &(E,) is the angular 
asymmetry: coefficient depending now on E,; both these functions are represented in 
Fig. 1. We see that æ varies in the limits from — 0.37 for slow neutrons to — 0.45 


"for fast neutrons, when for the case of a proton at rest in used approximation 


& = — 0.55. It means that the neutron asymmetry is considerably retained for 
muon capture in nuclei (67—80% of that in hydrogen). Then not only an experimental 
observation of this asymmetry may be expected, but the conserved vector current 
hypothesis may be probably checked too. 

5. Up to now there were reported three experimental results on neutron asymmetry 
in muon capture: 1) |gP,| = 0.03+0.15 [1], 2) « P, = — 0.05+0.02 [17] (both for 
capture in sulphur), and 3) œ — 4-0.154-0 18 [33] te capture in magnesium; this 
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Fig. 1. I (E,) — the neutron spectrum in arbitrary units, « (E,) — the angular asymmetry parameter. 


value was measured for emitted neutrons with E, > 5.5 MeV). These experiments 
are difficult because of the great gamma ray and evaporation neutron background, 
and their accuracy is for the present unsufficient to tell anything conclusive about 
the u--capture theory adopted here. 

The author wishes to express his cordial thanks to Dr. Chou Huang-Chao for 
collaboration and advice during his stay at the J.I.N.R.***) 
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